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Abstract 

We study the "hyperboloidal Cauchy problem" for linear and semi- 
linear wave equations on Minkowski space-time, with initial data in weighted 
Sobolev spaces allowing singular behaviour at the boundary, or with poly- 
homogeneous initial data. Specifically, we consider nonlinear symmetric 
hyperbolic systems of a form which includes scalar fields with a X(f> p non- 
linearity, as well as wave maps, with initial data given on a hyperboloid; 
several of the results proved apply to general space-times admitting confor- 
mal completions at null infinity, as well to a large class of equations with 
a similar non-linearity structure. We prove existence of solutions with 
controlled asymptotic behaviour, and asymptotic expansions for solutions 
when the initial data have such expansions. In particular we prove that 
polyhomogeneous initial data (satisfying compatibility conditions) lead to 
solutions which are polyhomogeneous at the conformal boundary of 
the Minkowski space-time. 
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1 Introduction 

Bondi et al. [6] together with Sachs [34] and Penrose [33], building upon the 
pioneering work of Trautman [36,37], have proposed in the sixties a set of 
boundary conditions appropriate for the gravitational field in the radiation 
regime. A somewhat simplified way of introducing the Bondi-Penrose (BP) 
conditions is to assume existence of "asymptotically Minkowskian coordinates" 
(x^) = (t, x, y, z) in which the space-time metric g takes the form 

„ _„ _ hju, (t - r, 6, ip) | hju, (t - r, 9, if) | 

&fj,u 'lnv — T r 2 "T ■ ■ ■ , 

where r]^ is the Minkowski metric diag(— 1, 1, 1, 1), u stands for t — r, with 
r,9,(p being the standard spherical coordinates on M 3 . The expansion above 
has to hold at, say, fixed u, with r tending to infinity. Existence of classes of 
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solutions of the vacuum Einstein equations satisfying the asymptotic conditions 
(1.1) follows from the work in [20] together with [3,4,18,19]. As of today it 
remains an open problem how general, within the class of radiating solutions 
of vacuum Einstein equations, are those solutions which display the behaviour 
(1.1). Indeed, the results in [1-4, 17] suggest strongly 1 that a more appropri- 
ate setup for such gravitational fields is that of poly homogeneous asymptotic 
expansions: 



In the context of expansions in terms of a radial coordinate r tending to infinity, 
the space of poly homogeneous functions is defined as the set of smooth functions 
which have an asymptotic expansion of the form 



for some sequences i%i,Ni, with rij f oo. Here the symbol ~ stands for "be- 
ing asymptotic to": if the right-hand-side is truncated at some finite i, the 
remainder term falls off appropriately faster. Further, the functions are 
supposed to be smooth, and the asymptotic expansions should be preserved 
under differentiation. 2 

The suggestion, that the expansions (1.2) are better suited for describing 
the gravitational field in the radiation regime than (1.1), arises from the fact 
that generic - in a well defined sense - initial data constructed in [1-4, 17] are 
polyhomogeneous. This leads naturally to the question, whether polyhomo- 
geneity of initial data is preserved under evolution dictated by wave equations. 
In this paper we answer in the affirmative this question for semi-linear wave 
equations, and for the wave map equation, on Minkowski space-time. We de- 
velop a functional framework appropriate for the analysis of such questions. 
We prove local in time existence of solutions for classes of equations that in- 
clude the semi-linear wave equations and the wave map equation on Minkowski 
space-time, with conormal and with polyhomogeneous initial data. We show 
that polyhomogeneity is preserved under evolution when appropriate (neces- 
sary) corner conditions are satisfied by the initial data. We note that the initial 
data considered here are more singular than allowed in the existing related re- 
sults [7, 28, 31]. We are planning to analyse the corresponding problems for the 
vacuum Einstein equation in a forthcoming publication, see also [30]. 

Our main results are the existence and polyhomogeneity of solutions with 
appropriate polyhomogeneous initial data for the nonlinear scalar wave equa- 
tion, and for the wave map equation. We achieve this in a few steps. First, 
we prove local existence of solutions of these equations in weighted Sobolev 

1 Cf. [29] and references therein for some further related results. 

2 The choice of the sequences n»,iVj is not arbitrary, and is dictated by the equations at 
hand. For example, the analysis of 3 + 1 dimensional Einstein equations in [17] suggests that 
consistent expansions can be obtained with m = i. On the other hand, Theorem 5.7 below 
gives actually n, = i/2 for wave-maps on 2 + 1 dimensional Minkowski space-time. We note 
that the 2 + 1 dimensional wave map equation is related to the vacuum Einstein equations 
with cylindrical symmetry (cf., e.g., [5, 14, 15]). 





(1.3) 



i=0 j=0 
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spaces, cf. Theorems 4.1 and 5.1. The next step is to obtain estimates on the 
time derivatives, cf. Theorems 4.4, 5.4 and 5.6. Those estimates are uniform 
in time in a neighbourhood of the initial data surface if the initial data sat- 
isfy compatibility conditions. Somewhat surprisingly, we show that all initial 
data in weighted Sobolev spaces, not necessarily satisfying the compatibility 
conditions, evolve in such a way that the compatibility conditions will hold on 
all later time slices; see Corollary 4.5 and Theorems 5.4 and 5.6. Finally, in 
Theorems 4.10 and 5.7 we prove polyhomogeneity of the solutions with poly- 
homogeneous initial data; this requires a hierarchy of compatibility conditions. 
We hope to be able to show in a near future that polyhomogeneity of solutions 
can be established, for polyhomogeneous initial data, with a finite number of 
compatibility conditions. 

The restriction to Minkowski space-time in Theorem 5.7 is not necessary, 
and is only made for simplicity of presentation of the results; the same remark 
applies to Theorem 4.1. Similarly the choice of the initial data hypersurface as 
the standard unit hyperboloid is not necessary. 

This work is organised as follows: First, the reader is referred to Appendix A 
for definitions, notations, and the functional spaces involved; we also develop 
calculus in those spaces there. In Section 2 we briefly recall Penrose's conformal 
completions, as they provide the link between the asymptotic behavior of fields 
and the local analysis carried on in this work. In Section 3 we consider linear 
equations. There the key elements of our analysis are: a) Proposition 3.1 and 
its variations, which give a priori estimates in weighted Sobolev spaces; b) the 
mechanism for proving polyhomogeneity, provided in the proof of Theorem 3.4. 
The transition from the linear weighted Sobolev estimates to their nonlinear 
counterparts is done in Sections 4 and 5. This has already been outlined above, 
and requires a considerable amount of work. In Appendix B we prove sev- 
eral auxiliary results on ODE's, some of which are fairly straightforward; as 
those results are used in the body of the paper in various, sometimes involved, 
iterative arguments, it seemed convenient to have precise statements at hand. 

Some of the results proved here have been announced in [16]. 

2 Conformal completions 

The aim of this section is to set-up the framework necessary for our considera- 
tions; the results here are well known to relativists, but perhaps less so to the 
PDE community. In any case they are needed to establish notation. Consider, 
thus, ann+1 dimensional space-time {J%,q) and let 



Let O h denote the wave operator associated with a Lorentzian metric h, 



~9 = tf Q . 



(2.1) 
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We recall that the scalar curvature R = R(q) of 3 is related to the corresponding 
scalar curvature R = R(q) of q by the formula 



m 2 = R-2n{ -O g n + __ L—M . (2.2) 




It then follows from (2.2) that we have the identity 

n+3 (_ „ n— l/^„2 



□0^-—/) = n-T- ^d b / + - R)f\ ■ (2.3) 

It has been observed by Penrose [33] that the Minkowski space-time (^#, 77) can 
be conformally completed to a space-time with boundary (^#,77), 77 = S7~ 2 77 
on by adding to ^ two null hyper surf aces, usually denoted by J? + and 
, which can be thought of as end points {^ + ) and initial points (<^~) 
of inextendible null geodesies [32,33,38]. We will only be interested in "the 
future null infinity" an explicit construction (of a subset of J^ + ) which is 
convenient for our purposes proceeds as follows: for (x ) 2 < >(£ 1 ) 2 we define 



M — 



(2.4) 



in the coordinate system {y^} the Minkowski metric 77 = — (dx ) 2 + (dx 1 ) 2 + 
(dx 2 ) 2 + (dx 3 ) 2 = T] a pdx a dx 13 takes the form 



r? = ^ a pdy a dyP , n = 77^7/V . (2.5) 



We note that under (2.4) the exterior of the light cone Cq M = {r} a px a x^ = 0} 
emanating from the origin of the x^-coordinates is mapped to the exterior of 
the light cone Cq = {Tj a py a y^ = 0} emanating from the origin of the y^- 
coordinates. The conformal completion is obtained by adding Cq to , 

with the obvious differential structure arising from the coordinate system y^ . 
We shall use the symbol J to denote Cq \ {0}, and to denote Cq \ {0} n 
{y° > 0}. As already mentioned, so defined is actually a subset of the usual 
r y, but this will be irrelevant for our purposes. 

We note that (2.4) is singular at the light cone Cq M . This is again irrelevant 
from our point of view because we are only interested in the behaviour of the 
solutions near J? + , and finite speed of propagation allows us, for that purpose, 
to disregard what happens near Cq M . 

The above procedure can be adapted for several metrics of interest, such 
as the Schwarzschild, Kerr, or Robinson- Trautman metrics, to similarly yield 
conformal completions of space-time by the addition of null hypersurfaces 
This observation was at the origin of Penrose's proposal to describe systems 
which are asymptotically flat in lightlike directions through the use of conformal 
completions. 
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It is noteworthy that the conformal technique allows one to reduce global- 
in-time existence problems to local ones; this has been exploited by various 
authors [8-13] for wave equations on a fixed background space-time. Further, 
Friedrich [22, 23, 26] has used this approach to obtain a global existence result 
for Einstein equations to the future of a "hyperboloidal" Cauchy surface, with 
"small" smoothly conformally compactifiable initial data, cf. also [21,24,25]. 

On a more modest level, the identity (2.3) can be used as a starting point 
for the analysis of the asymptotic behaviour of solutions of the scalar wave 
equation near J? + , as it reduces the problem to a study of solutions near a null 
hypersurface. This is the approach used in this paper. There are associated 
identities for fields of any spin [33], which provide a convenient framework for 
similar questions for those fields. 

3 A class of linear symmetric hyperbolic systems 

In this section we define a class of linear symmetric hyperbolic first order sys- 
tems on a set of the form M XQ x /, where / is an interval corresponding to 
the time variable, which will be denoted by r, and we derive our key energy 
inequality in weighted Sobolev spaces. (We note that in some of our further 
applications the vector d/dr will be lightlike, and not timelike as is usually 
the case. It should be pointed out that in our conventions the time variable 
is the last coordinate, allowing x to be the first variable, consistently with the 
conventions of the preceding sections.) We start by introducing some notation 
for the sets within the "space-time" M XQ x /, which will be relevant in what 
follows 3 : 

t > 0, 2(x 2 + 1) < x 1 < x , ^x 2 ,x!,t = {r = t, x 2 < x < xi - 2t}, 

(3.1a) 

T > 0, 2(x 2 + T) < Xi < X , ttx2,xi,T = Uo<r<T ^x 2 ,x u r, (3.1b) 

< 2t < xi < x , S xijt = {r = t,Q < x < x ± - 2t}, (3.1c) 
< 2T < xi, Q xuT = Uo<t<T s zi,t • (3.1d) 

There is a natural identification between £ X 2,xi,t and M X2X1 —2t, similarly be- 
tween T, Xlt t and M El _2t, and we shall freely make use of such identifications 
throughout. We shall write ||/(t)||.#»« for \\fh X2 , xl J\^(i: X2 , xl , t ), or for ||/ Is^, J^s^), 
etc.; the distinction should be clear from the context. 

3 The motivation for the factors of 2, and the general form of the sets considered, arises as 
follows: The set dM x / should be thought of as a smooth null hypersurface in space-time; 
e.g., in Minkowski space-time with Minkowskian coordinates y M , it can be the intersection 
of the half-space {y° > i} with the light cone emanating from the origin — . Then 
t is the Minkowski time, perhaps shifted by a constant, say r = y° — |. The coordinate x 
is a coordinate which vanishes on dM x I, in the current example e.g. x = \J^2(y 1 ) 2 — y°- 
Finally, in such a Minkowskian setup, the hypersurfaces x = xi — 2r, which determine one of 
the boundaries of the S's and fi's denned in (3.1), correspond to the converging light cones 
V° + \f~Yl(y 1 ) 2 ~ const. The restrictions 2(2:2 + t) < xi < xo (in the definition of T, X2iXli t) 
and 2(x 2 +T) < x\ (in the definition of CI X2 ,x 1 ,t) are not necessary, and are only made for 
simplicity of discussion. 
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We shall be interested in symmetric hyperbolic first order systems which in 
local coordinates take the form 



[A^)d li + A{z)]f = F, (3.2) 

where z v = (y l , r), with the following properties: 

^1) / and F are sections of a bundle which is a direct sum of two N\ and 
N2 dimensional Riemannian bundles; we will write 

'=(;)• f -{:)- ^ 

In local coordinates ip and a are thus W Nl valued, while ip and b are R^ 2 valued. 
The respective scalar products will be denoted by (•, -\ and (-, -) 2 . We shall use 
the generic symbol V to denote 4 a covariant derivative compatible with those 
scalar products, e.g., if X is a vector field on SI XOj t, then 

X ({<(>, ^ ) = (Vx0, + {</>, Vx ^ , (3.4) 

similarly for (-, -) 2 . V will also be assumed to be compatible with every other 
structure at hand whenever useful in the context, e.g. a Riemannian metric on 
M, etc. 

^2) The left hand side of (3.2) can be written as 

ti>_V, V +14 \ + (Bn ,3. 5) 



where L is a first order differential operator. Here denotes the formal adjoint 
of L, in the sense that if f2 = M, or M Xl , or M X2)X1 , and if tp,ip are in Ci(J7), 
then 

/ (<P,Lip\ dfi= (tfip,ip} 2 dfi , (3.6) 
Jn 

where dfj, is a measure on M which will, we hope, be obvious from the context. 
By density Equation (3.6) will still hold with f2 = M X2tX1 for all a,f3 £ R, 
all if £ J^ a (M X2)a;i ) and all ^ £ ^1 (-^X2,xi)- Equation (3.6) forces L not 
to contain any r- or 2- derivatives, where the letter x denotes a coordinate as 
defined in Section A, thus 

L = £ A (x,v,T)d A +£(x,v,T) . (3.7) 

It follows that the principal part of the system (3.5) is of the form 

/ E»_d» i A d A \ 

\ (£ A Yd A E»d» ) ' [6 * } 

where A t denotes the transpose of a matrix A. Equation (3.8) explicitly shows 
that (3.5) is symmetric hyperbolic when the £±'s are symmetric with E± pos- 
itive definite; the notions of "symmetric hyperbolic" and "symmetrizable hy- 
perbolic" are identified throughout this work. 

The hypotheses above will be assumed throughout this section. 

4 In some situations (3.4) might fail to hold, and some undifferentiated supplementary terms 
will occur at the right-hand-side of (3.4). We note that our results will not be affected by the 
occurrence of such terms, provided those terms satisfy bounds as in (3.17). 
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3.1 Estimates on the space derivatives of the solutions, a < 

-1/2 

Let us pass now to the description of the hypotheses needed to derive weighted 
energy estimates for space derivatives of /. To obtain such estimates, we shall 
require the existence of a constant C\ such that the (matrix-valued) coefficients 
i A and £ satisfy, in the relevant range of r's, 

ll^)ll Wi _ 2t) + E \\t A (r)yo {Mxi _ 2T) < a . (3.9) 

A 

Similarly writing 

Lt = tfAfa ^ t ^q a + £ t( x> Wj r ) ; ( 3 10 ) 

we require 

¥\r)yo {Mxi _ 2T) + E II^V)II Wi _ 2t ) < d . (3.11) 

^3) The matrices E± are symmetric and satisfy 

£±n M > eld , E^dpx < -eld , \E^x\ < C lX , (3.12) 

for some e > 0. Here n M denotes the field of future directed (i.e., 6(dr, n) > 0) 
6- unit normals to the surfaces {r = const}, where b is an auxiliary Riemannian 
metric h on M. (Later on we will mainly be interested in the case of E+s of 
the form E^_ = ® Id, for some vector fields e±.) For simplicity we shall also 
assume 

diE T ± = ; (3.13) 

this is by no means necessary, but is sufficient for the purposes of this paper. 
We will further assume 5 that the Et 's satisfy a bound of the form: 

(3.14) 

As far as the E^s are concerned, we allow singular behaviour which should, 
however, be somewhat less singular than 1/x; to control that, we require ex- 
istence of a function ( : R + — > R + , satisfying lim^o C( x ) = 0, such that for 
< x < X\ — 2r we have 

II^(t)||^-i (Mx) + \\El(r)\yo {Mx) + ||sV^(r)|| L =o (Mie) < ((x) • (3.15) 
When the operators E±V^ are written out explicitly as 

E&v = + B± , (3.16) 

5 We use a convention in which the covariant derivatives V M i5^ include terms associated 
with the vector density character of X M defined by (3.21); in particular this should be taken 
into account when verifying that the estimates (3.14)-(3.15) hold. 
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we require that 



\\B-{t)\\^ [Mxi _ 2t) < d , ||B+(T)||^-i (AfB) <C(a;), 0<x< Xi -2t. 

(3.17) 

^4) The matrices B a f,, a,b = 1,2, satisfy the bounds 

\\ B ^ t )W^\m x1 . 2t ) + \\ B ^ t )\\^'\m x1 _ 2t ) + ll^nWII^^-^) < Ci , 

\\B 22 {r)\\^ (Mx) <ax), (3.18) 

this last equation holding again for < x < x\ — It. 

Our final hypothesis concerns the "acausal" nature of the boundary of 

x 2 ,x!,T is "non-timelike", in the sense that for any covector n^, which 
is positive on outwards-pointing vectors and vanishes on vectors tangent to 
dttx2,x u T we have, on dQ X2jXljT n {r > 0}, 

££n„ > . (3.19) 

(We note that (3.12) already guarantees that (3.19) holds on dQ X2tXl ,T n {r = 
T or x = 0}.) 

The essential point of the above hypotheses is that the boundary {x = 0} 
is characteristic for Equation (3.2), with the dimension of the relevant kernel 
being constant over the boundary. 6 Weighted estimates, in the spirit of Propo- 
sition 3.1 below, near such characteristic boundaries hold for general symmetric 
hyperbolic systems, this will be discussed elsewhere. 

Weighted energy inequalities in Jif^ spaces with arbitrary values of k may 
be proved under various hypotheses on the coefficients which appear in (3.2). 
We note one such result for systems satisfying ^l)-^5), which lies in line with 
our remaining investigations. The restriction a < —1/2 seems to be inherent to 
the problem at hand. We consider first the case a < —1/2; the case a = —1/2 
is handled by the same methods, under somewhat more restrictive conditions 
on the coefficients, in Section 3.2. 

Proposition 3.1 Suppose that a < —\, k > ^ + 1, k G N, and set either 
/(*) = f\v xl , t , < X! < xq, < t < t max = xi/2, or f(t) = /|e X2 ,^, 
< 2x2 < x\ < xq, < t < t max = x\ — 2x2- Under the hypotheses *^l)-^5), 
there exists a constant C2 depending upon x\, C\, n, N, k and a, as well as 
upon the "error function" ( and the boundary manifold dM, such that for all / 
satisfying (3.2) for which /(0) G H l ° c and for all < t < t max we have 

with an identical estimate with M X1 _* replaced by M X2)X1 _*. 

6 We are grateful to H. Friedrich for useful discussions concerning this point. 
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Remark: The condition fc>n/2 + lis needed to obtain Ci-weighted control 
of the solution; there are no restrictions on k if we have at our disposal an 
a priori C\ weighted bound for /, and if the coefficients in the equation are 
suitably regular. In such a case, for k < n/2 + 1, the inequality (3.20) should 
be modified by adding a term ||/(s)||^a( M 2 ) under the integral appearing 
in (3.20). 

Proof: We start by proving (3.20) on sets M X2jXl - t ; in that case we are mainly 
interested to obtain uniform control for small values of X2, with eventually X2 
tending to zero; without the uniformity in X2 the estimate would of course 
be standard. Keeping this in mind, let be the "energy-momentum vector 
density" , 

X»= x- 2a - x ^{^^>,E^^\ + ^^,E%^i)\}. (3.21) 

0<|/3|<fc 

Suppose, first, that /(0) G H l £?^ standard results [35, Vol. Ill] show that 
f(t) 6 -fffc+i, and we then have 5 

V M X" = N 1 +D 1 +D 2 + E 1 + E 2 + E 3 , (3.22) 

where 

0<\/3\<k 
0<\/3\<k 



0<\8\<k 



El = (2A -2a- l)x- 2a - 1+2 ^(^<p , {E - d » x) &y\ , 

0<\f3\<k 
0<\(3\<k 

E 3 = x- 2a - 1+2/3l (^,(V M ^)^) 2 • (3-23) 

0<\(3\<k 

Since 2a + 1 < 0, from (3.12) one finds that 



It,-. 



Nidxdv < -\2a + l\e\\ipf ,i (3.24) 



which is strictly negative except if is identically zero, and can be used to 
control some of the error terms which occur at the right hand side of (3.22). 
(Here we have used the form (A. 4) of ||Y>|| 2 . .) For example, to control E% 



we take any X3 satisfying 2x2 < ^3 < x\ — 2s (we will make a more precise 
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choice of X3 later), and we write 

/ E ?> dxdv = £3,1 + £3,2 



£3,1 = / E 3 dxdv , 

J^x 2 ,x 1 ,s C\{ x > x a} 

E 3: 2 = E 3 dxdv . 



By (3.15), £3,2 can be estimated as follows: 

\E 3 , 2 \ < Yj I C{x)x- 2a - 2+2 ^\^ij\ 2 dxdv 



0</3<k x 2' x i' a C\{ x < x 3} 

< ^±^H\\ 2 



10 

if X3 is chosen small enough. Once this choice has been done, we can clearly 
estimate E 3 1 as 

£3,1 < c\m\ a , 

with some constant which is determined by x 3 . The integrals of the error terms 
Ei and E 2 are estimated in the obvious way, cf. (3.12) and (3.14): 



J {E x + E 2 ) dx dv < C\\ip( 

J ^Xo .11 .s 



k 



- J X2 , s 

To control the terms Di and D 2 we use the evolution equations (3.5): 
E^V^if = ^(^V^>) + [£/%,i^> 

= -L^ip + ®P a + E{ , (3.25) 
E% = L]ip + [EtVp, 2% - @P{B lW + B 12 il>) , 

E^V^il) = L^^ip + ^b + E^, (3.26) 
El = tf\<p + [^V M , - ®\B 2W + B 22 ^) . 

Integrating D\ + D 2 over S :E2iXljS , one finds that the terms containing L&^ip 
and —L^St^ip in (3.25) and (3.26) cancel out; the terms containing S)^a and 
are estimated as (here the somewhat arbitrarily chosen factor 10 can be 
replaced by any other larger number if desired) 

2 ^ / 3.-20,-1+201 ^gP^gPa), + (^V,^> 2 ) 

dv 

<l|y||^ + ll«ll^+ (2a + 1)£ |I^H 2 + i + 7 ll&ll 2 1. 

- \vr\\3ie k 11 11^ 1Q imi^«+i (2a + i) e n n ^-i 
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The terms containing the commutators [@P,L]ip and [@P,L^\ip, can be esti- 
mated 7 using the weighted commutator inequality (A. 35), while the Bn, -B12, 
etc., terms can be estimated using (A. 34), by an expression of the form 

CCr (\\nle« + IMl3r- + (2 ^V )£ |MI^ + i) • (3-27) 
To estimate the commutator terms arising from E±, we calculate, e.g. 

k 

x k [E^d„d k } X = Y.^V{diE^)x k - i d k x -%x 
i=i 

= Eq + e 7 , 

k 

The terms arising from Eq are estimated in a straightforward way as in (3.27) 
using (A. 36). The dangerous term £7 can be written as 

k k 

E7 = £(|K(4£±)* fe -^ +1 X = Y,&)x i -\di- 1 d*El)x k - i+1 d k - i+1 X , 
i=i i=i 

and can thus again be estimated as in (3.27) provided that d x E x _ G Sf^, that 
d x E+ G ^T 1 !, and that (3.15) holds. Other terms in the i£± commutators are 
handled in a similar way. 

Summarizing, we have derived 



L 



< Cd \\a( S )\\ 2 ^ + \\b(s)f a _ h + ||^( S )||3ej- + ||v( S )||5r- 



(3.28) 

where d n // stands for dxdv, or for any measure uniformly equivalent to dxdu. 
Stokes theorem, 

/ V^dydr = j X^dS^, 
and our hypotheses on the geometry of the problem lead to 

||/(t)||^« < C (||/(0)||^ + d £ (||a(s)||^« + ||6(s)||^ rl/2 + ||/( S )||^) ds 

Gronwall's lemma establishes (3.20) on the family of hypersurfaces £ X2 ,xi,t for 
/(<) G tf^. If /(t) G tf£ oc , we approximate /(0) by a sequence of functions 
/ n (0), with / n (0) G converging to /(0) in ^. a (Sx 2 ,xi,t), and we solve 

Equation (3.2) with initial data / n (0). The inequality (3.20) applied to the 

7 This step requires weighted L°° control of <j> and tp, and weighted W 1,oa control of the 
coefficients in the equation. The hypothesis k > n/2 + 1 is not needed if such a priori bounds 
are known. 
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functions f n (t) — f m (t) shows that f n (t) is Cauchy in passing to the limit 
n — ► oo the desired result for /'s such that /(0) G J^, Q (S :E2)Xl) j) easily follows. 

Since all the constants above are xi independent, an elementary argument 
using the monotone convergence theorem shows that the inequality (3.20) for 
the S X i,t's follows from the one for the S X2jXljt 's by passing to the limit X2 — > 0. 
□ 



3.2 Estimates on the space derivatives of the solutions, a = 

-1/2 

When a = —1/2 we do not have the f3± = negative terms in Aq at our dis- 
posal in Equation (3.23), so that we cannot allow coefficients as singular as in 
the previous section. To handle that case we keep all the structure and regu- 
larity conditions already made, with the following supplementary restrictions: 
Equation (3.15) is replaced by 

\\Ei(r)\yo {Mx) + \\E x + {r)\\ n{Mx) + || V^(r)|| L =o (Jl#ie) < d . (3.29) 

Instead of (3.17) we require that 

||£±(t)||^o (M:ci _ 2t) <C' 1 , (3.30) 
while condition (3.18) becomes 

WBabirn^^KC^ (3.31) 

We then obtain: 

Proposition 3.2 Suppose that fc > f + 1, fc G N, and set either f(t) = f\n x t , 
< xi < x , < t < i max = xi/2, or f(t) = /|s X2 , xl , t , < 2x 2 < x x < x , 
< t < t max = x 1 - 2x 2 - Under the hypotheses ^l)-^) together with (3.29)- 
(3.31) there exists a constant C2 depending upon x\, C\, n, N, k and the boundary 
manifold dM, such that for all / satisfying /(0) G H l k oc and for all < t < t max 
we have 



e c ^- s H\\a(s^ 2 



(M Xl -2s) 



"wiiW^-) - ° 2eCs, ( ll/(0)ll V" 2 («.,) + i 

+ ■ (3 - 32) 

with an identical estimate with M X1 _* replaced by M X2 ^ X1 -.*. 

Proof: The proof is essentially identical, but simpler, to that of Proposi- 
tion 3.1. We simply note that the key inequality (3.28) gets replaced by 



< Cd (\\a(s)\\^_ 1/2 + \\b(s)\\^_ 1/2 + ||^)||^- 1/2 + IbtollJ 



(3.33) 
□ 
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3.3 Estimates on the time derivatives of the solutions 



The hypotheses assumed in the previous section ensure that we can algebraically 
solve Equation (3.2) for d T f, and then recursively obtain formulae for d l T f . 
Under the hypotheses of Proposition 3.1, it is then straightforward to obtain 
estimates on the norms 

\\((xd T yf)(T)\u ka _^ xi _ 2T) , o< i <k, 

provided suitable weighted conditions are imposed on the r derivatives of the 
coefficients of Equation (3.2). However, we would like to obtain derivative 
estimates without the x factors, uniformly in r. Clearly a necessary condition 
for the existence of such estimates is that 

ll(4/)(0)||^_ i (E, 1 )<oo, 0<i<k. (3.34) 

It turns out that (3.34) does not need to hold for arbitrary initial data /(0) £ 
J^l? 1 , and the requirement that it does lead to the j-th order compatibility con- 
ditions: by definition, these are the conditions on f(0) which ensure that Equa- 
tion (3.34) holds for < i < j. Since, for sufficiently differentiable solutions 
of Equation (3.2), all the derivatives d T f(0) can be explicitly written as an i- 
th order differential operator acting on /(0), the compatibility conditions are 
conditions on the behaviour of the initial data /(0) near the "corner" x = 0; 
we shall therefore sometimes refer to them as 11 corner conditions" . We note 
that there can be corner conditions in weighted Sobolev spaces, or in weighted 
Holder spaces; in this section we will be mainly interested in the latter, defined 
by Equation (3.38) below. 

The following example is instructive in this context: For < t < y let g be 
a solution of the 1 + 1 dimensional wave equation 

d 2 d 2 \ 

SS- ^5)9 = 0, (3-35) 



with initial condition 
9 



dt 2 dy 2 



tt+i d 9 



= 2cy . 

t=o at 



t=o = 2(a + l)y a , 



for some constants C,a£l. From Equation (3.35) we can obtain a system of 
the form (3.5) by introducing r = t, x = y — t, ip = (g, (d T — 2d x )g), if) = d T g, 
and setting L = 0, E^d^ = d T ® id, E+d^ = (d T — 2d x ), so that we have 



3r 



9 \_(*\ = ( G 

(d T - 28 x )g ) V ) V o 

(d T -2d x )1> = 0. 



The solution is 



= ( C + l)(y + tr +1 + (C-l)(y-tr +1 
= (C + 1)(2t + x) a+1 + (C - l)x a+1 . 



14 



It follows that for each < r < 1, k G N, and [3 < min{0,a + 1}, we have 
g(r, •) G J4?f((0, 10]), consistently with Proposition 3.1. Somewhat surprisingly, 
for r > and for all « G N the functions <9* <?(r, •) are smooth in x up to x = 0. 
However, the L°° bound for <9*g(r, •) blows up as r tends to zero except in the 
case 

C = -l. (3.36) 

Condition (3.36) is precisely the corner condition needed for d T g(0, •) to be 
better behaved than d x g(0, •) at x = 0. In the example under consideration the 
fulfillment of the first order corner condition guarantees already that all the r 
derivatives of g will be well behaved, but we do not expect this to be true in 
general. 

Let us pass to a derivation of the desired estimates. We shall use a method 
which works directly in weighted Holder spaces, avoiding the use of weighted 
Sobolev spaces; the price one pays is the need to consider systems somewhat less 
general than (3.5), but still general enough for our purposes. More precisely, in 
this section we restrict our attention to systems of the form 

d T <p + Buip + Bi 2 i> = L n (p + Li2^ + a, (3.37a) 

e+Y> + B 2 i<p + B 2 2^ = L 2 i<p + L 2 2tp + b , (3.37b) 

with 

e + tf) = (d T - 2d x )i) . 

We assume that the L a bS, a,b = 1,2 are first order differential operators of the 
form 

L ab = L^d A + xU ah 3 T + xL%d x , (3.38) 

with bounded coefficients L^ b ; no symmetry hypotheses are made. Clearly the 
intersection of systems of equations satisfying (3.37) with those of the form 
(3.5) is non-empty. (As we will see in Sections 4 and 5 below, non-linear wave 
equations on Minkowski space-time can be written in the form (3.37).) In 
particular Proposition 3.1 provides a large class of solutions of (3.37) such that 

(<p,ip){r) G J? k a {M xl _ 2T ) C % a (M xl _ 2T ) 

for £ < k — n/2. We shall therefore assume that a solution / = (ip, ip) satisfying 
/(t) G ^ a (M Xl _2r) is given, and study its r-differentiability properties. For 
the purposes of the proof below it is convenient to introduce auxiliary spaces 
^ p (Cl) defined, for p < £, as the space of functions / in Cg(fl) such that the 
norm 

ll/Kyn) = sup x- a \{xd x )\xd T y^ v d k T f\ 
0<i + j + k + \j\ <l 
0<k<p 

is finite. Obviously, ^ £ = ff". Similarly one defines «jf (0) using the norm 

H/ll^,, = sup (i+\^x\)^x- a \(xd x y(xd T y^d k T f\ . 

0<i + j + k + \j\<£ 
< k < p 
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Clearly ^(fi) = We shall write ^ for <^/. 

Proposition 3.3 Let a < 0, I G N, write f2 for Q X1j t (with Q X1j t as in (3. Id)), 
and suppose that L^ b ,B ab G <^°(fi), a G ^(^j, & G ^-h n )- Consider 
f = (ip,ijj) — a solution of (3.37) satisfying 

Vre[0,T] /(r) G <^(M xl _ 2r ) . 

Then: 

1. For all e > we have 

M)€^ e f 2i (nn{x + 2 T >e}) . 

Further, for any r > the compatibility conditions of order p = [£/2\ (the 
integer part of £/2) are satisfied by (</?(t), V>( T )) : 

Vl<Kp ^¥>(r), 4V(r) G (M Xl ) , (3.39) 

Here /3 = L^/ 2 J if a = 0, and /? = otherwise. 

2. If there exists 1 < p < £/2, p G N, such that Equation (3.39) holds with 
/3 = at r = 0, then 

(¥>^)e^4(n)c^(n), (3.40) 

with /3 = p if a = 0, and = otherwise. 

Remark: The method of proof here gives a number of well controlled time 
derivatives smaller by a factor 2 than the number of space ones. This is, how- 
ever, irrelevant, when £ = oo, which is the main point of interest in this work. 
We note that energy estimates as in the proof of Theorem 4.4 below provide an 
alternative, more complicated way of establishing a stronger statement, with 
more controlled time derivatives for large ts. In the linear case considered 
here the function F occuring there vanishes, so that all the complications aris- 
ing from the non-linearities disappear, and somewhat stronger results can be 
obtained using the methods there. 

Proof: By rearranging terms and redefining the L^'s, the B^s, and the source 
functions a and b we may without loss of generality assume that 

Kb = o • 

One can rewrite Equations (3.37) as xd T ((p,ip) = a partial differential oper- 
ator linear in xd x and d v ; by iteration this immediately yields (<p, ijj) G 
Equation (3.37a) shows then that d T ip G ^_^ , hence tp G ^ r On the other 
hand, Equation (3.37b) gives e+(ip) G ^L^o + hence d T e + (ip) G «^f . 

In order to extract further information out of Equation (3.37b) we use the 
straightforward identity 

pt+x/2 

^(x,v a ,t) = tp(x + 2t,v a ,0) + / e + (ip)(2v,v A ,r -v + x/2) dv . (3.41) 

Jx/2 
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(We note that for each e > the first term above is uniformly Q on the set 
n {x + 2r > e} n {x < xq}.) Differentiating Equation (3.41) one obtains 



rr+x/2 

d T ^(x,v A ,T)=d T ip(x + 2T,v A ,0)+ / d T e + (i;)(2v,v A ,T-v + x/2) dv 

Jx/2 



since a < and d T e + (ip) G ^^jo' straightforward estimations show that d T ip G 
^-2|o> hence ^ G if a / 0, while ^ G when a = 0. 

Let /3 r = if a / and j3 r = r when a = 0, and suppose that 99 G ^f{_ T ^ T 

and V G «£jf r for some 1 <r < (I — l)/2; we have already shown this to hold 
for r = 1. Equation (3.37a) gives 

It then follows from Equation (3.37b) that 

e+W€^Li|r =► ^M^C-V 
Differentiating r + 1 times Equation (3.41) with respect to r we obtain 

rr+x/2 

d r T +1 ^(x,v A ,T) = d r T +1 ^(x+2r,v A ,0)+ / d r T +1 e+{ip){2v, v A , t-v+x/2) dv , 

Jx/2 

which gives d r T +l ip G ^^2r-2|0' nence ip G ^g^r-Hr+v anc ^ tn e induction is 
completed. □ 



3.4 Poly homogeneous solutions 

We now wish to show that solutions with poly homogeneous initial data will be 
polyhomogeneous. Let Q X0 ,T be defined by Equation (3. Id); we shall denote by 
£/£(£1 XOj t) the space of functions / defined on Q Xo ,T which can be written in 
the form 

k Ni 

^2 ^2 X l& \n 3 Xfij + f a+ kS+e , 
i=0 j=0 

for some e > 0, some functions fij G C 00 (Q XQj t), and some sequence (JVj) of 
non-negative integers. We also require that f a+ kS+e £ < ^'So' kS+€ {^xo,T)- We set 

:= n fceN ^f . 

The following properties are useful in what follows: 

• If < x\ < xq — T/2, then a function / G C oo (Q X0! t) is in ^(VL Xo t) if and 
only if for any coordinate patch C of dM we have / G &/ k s (U xl ), where 
U XI =]0,xi[x^ x [0,T], and if / G Coo(^t), where O int = O^r n {x > 

Xi}. 

• For all e > we have tf£> +pS+e C xP*/*\ in particular ^ C 
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• It does not hold that g/£ C c &® , however, for all e > we have srf k 5 C 
c €^. More precisely, if / € then there exists p G N such that 
{l + \lnx\ 2 )-P/ 2 f G<T£. 

• As before we assume that 1/5 G N, which implies £=c^ C ^+i/S c ^fc+i ■ 

• is stable under multiplication: if G then /</ G . 

• ^ is stable under differentiation with respect to r and to u, as well as 
under xd x : if / G *f£, then <9 T / • / (i > 2), zd*/ G ^f, with the 
vector fields X; L defined in Section A, cf. Equation (A. 7). 

In this section we will consider systems of the form 

d T ip + Bnip + B 12 tp = L U ip + L 12 i/j + a , (3.42a) 
d x tp + B 2 i<p + B 22 ip = L 21 ip + L 22 ?p + b , (3.42b) 

with the Ly's, i,j = 1,2 of the form 

Lij = Lp A + Ljjdr + xL?jd x , (3.43) 

with 

^exV^, L^,L^,L^ 2 e^ . (3.44) 

No symmetry hypotheses are made on the matrices L^-. Conditions (3.42a)- 
(3.44) are easily seen to be compatible with those made elsewhere in this paper, 
cf., e.g., the proof of Corollary 3.5 below. The reader is warned, however, that 
the operators here do not coincide with those in (3.37): to bring (3.37) into 
the form (3.42) one needs to multiply Equation (3.37b) by —1/2, transfer the 
operator d T from the left- to the right-hand-side of (3.37), and appropriately 
redefine the L 2 j : s. 

We start with the following result, which assumes that the solutions have 
both space and time derivatives controlled, in the sense of weighted Sobolev 
spaces; recall that this hypothesis can be justified for equations satisfying more- 
over the hypotheses of the previous sections: 

Theorem 3.4 Let (3, (3' G E, k G NU{oo}, and let ((p,ip) be a solution of (3.42) 
in ^)(^xo,t)- Suppose that (3.44) holds, and that 

B u G (*/* n L°°) (fi X0)T ) , B 12 , B 22 , B 21 G <(fi* 0) T) , (3.45a) 
a, b G x^ k \n xo , T ) , <p(0) e x^ k 5 (M X0 ) . (3.45b) 

Then 

cp G (x P ^i + ^5f) (^ ,t) , Y> G (^x ,t) + Coo(^^ ,t) • 

If one further assumes 

Lf 2 ,B 12) a ) #er(fi IOiT ) , 

then it also holds that 

v?g (x^ + 4 j nL°°) (n XOiT ). 
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Proof: It is convenient to decompose fin in the obvious way as 

fin = fi?i + fin , 
with B s u G x b ^_ x and fi n G C^. We rewrite (3.42) as 

d T ip + B^ip = a , (3.46a) 

d x ^ = ca, (3.46b) 

where 

c\ := Lu<p + L 12 ip + a - Bi 2 tp — B^ip , (3.47a) 

c 2 := ^21^ + -^22^ + b - B 2 if - B 22 tp , (3.47b) 

In what follows we let e > be a positive constant, which can be made as 

small as desired, and which may change from line to line. We note that C2 is in 
c €^o~ <L + x 13 ^, and integration in x of (3.46b), together with Propositions B.3 
and B.6, gives 

1p = Y> + Y>/3'+l-e + V'phg , 



where 

V'o(-) 



lim x _o ip(x, •) , if ft' + 1 - e > 0, 
, otherwise, 

with 
hence 

^ G Goo + + s0+ V fc 5 . 

Since Lj^ G ^' +<5 ~ e (d^ G < r£' 1 and xL\ x G x^ 5 n V$ C similarly for 
the other derivatives), we find that 

ci G + x'V/ + . 

We can then apply Proposition B.4 to (3.46a) to conclude that 

<p (=*/* + + ^ +p5 ~* , (3.48) 

with p = 1. Coming back to c 2 we find now that c 2 G + x 13 ^ + 't?£ > +p5 ~ e , 
and by Proposition B.6 we obtain 

iji G Coo + + x^+V/ + ^' +p,5+1 - e , (3.49) 

still with p = 1. To conclude, we proceed by induction; let ft' + p<5 < /3 + k and 
suppose that Equations (3.48)-(3.49) hold; it follows that c x G + x 13 ^ + 

< ^> +< ' P+1 ' ) ' 5 e • Applying Proposition B.4 to (3.46a) gives (3.48) with p replaced 
byp+1. It follows that c 2 G ^ +x /3 +^ +{p+1)5 ~ t ; Proposition B.6 applied 
to (3.42b) gives (3.49) with p replaced by p + 1, and the result is established. 
□ 

As a straightforward consequence of Theorem 3.4 we obtain: 
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Corollary 3.5 Let (3' G R, let (cp,ip) G ^'(^.t) be a solution of the system 
(3.5), and suppose that 

B ij ,E£,B±,£,tf,£ A ,(£ A y G ^ k S (n x0:T ) , (3.50a) 
Eland El — invertible, with (EL) -1 , (E'D^ 1 G ^(^I xo ,t) , (3.50b) 
(EL)- 1 ^ G x « n ^) (n XfhT ) , (el)- 1 e a G ^^f_ 1 (fi X0)T ) , 

(3.50c) 

(ELrHBn + BJ) G L°°(n X(hT ) . (3.50d) 

If 

a, b G s^(n X0)T ) , ^(0) G x^(M, ) , 

with /? G R, then 

^ G (xPj/s + ^ ( ^ ot) ) ^ G (>+ V* + x ^f) (n X0)T ) + Coo(n^) . 

In particular, if fe = cx> then the solution is polyhomogeneous. 

Proof: : We write Equation (3.5) as 

d T ^ + (E^y 1 {(Bn + b_)p + #} = (^D-^^i^-^a^ + a) 
- l^)" 1 {^V - (^22 + 5+)^} = (ED^a^^ + ^^ + ^AV + ft), 

(3.51) 

which is of the form (3.42), and we note that the hypotheses made on the 
coefficients of Equation (3.51) imply those of Theorem 3.4. □. 

An unsatisfactory feature of results such as Theorem 3.4 is that uniform 
estimates both on space and time derivatives of the solutions are assumed. 
Recall that uniform control of time derivatives can be obtained only if corner 
conditions are satisfied, and the hypotheses of Theorem 3.4 require an infinite 
number of those to be fulfilled. The same techniques can be used to obtain 
various expansions of solutions when a finite number of time derivatives are 
controlled only, but the statements turn to be out somewhat less elegant. We 
give an example of such results when 5 = 1: 

Theorem 3.6 Let (3 G R, k G N U {oo}, and let (ip,tp) be a solution of (3.42) 
in #f (O X0)T ) for some £ > 1. If Equations (3.44)-(3.50) hold with 6 = 1, then for 
any A < 1 we have 

v g (VV, 1 + + n,_ 2j _ 2 > <:g+ 2 ) (JW) > 
V* g (V+Vfc 1 + x^l + n<_ 2j _i> tf£lj+i +A ) (n X(hT ) + CooOQ (3.52) 

If one further assumes 

L» 2 ,B 12 ,a,<p(0)eL°°(n xo , T ) , 

then it also holds that 

tp g (x'H 1 + < n l°° + n^_ 2j _ 2 >o<_ + 2 ^2) (n X0)T ) • 
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Proof: The result is obtained through a repetition of the proof of Theorem 3.4, 
keeping track of the differentiability of the remainder terms. □ 

We are ready now to prove polyhomogeneity of solutions of the Cauchy 
problem for Equation (3.5). We consider only the simplest case of equations 
satisfying the conditions (3.53) below, considerably more general statements 
can be proved using similar methods. The differentiability hypotheses below 
are clearly satisfied by equations with smooth bounded coefficients; however, 
they also allow for a wide class of equations with poly homogeneous coefficients. 
We restrict ourselves to the case in which the corner conditions are satisfied to 
arbitrary order; if not, one obtains expansions as in (3.52), with a remainder in 
which a finite number only of time derivative are controlled; such results can be 
proved by identical arguments, compare the proof of Theorem 3.6. We hope to 
be able to show in a near future that the corner conditions are not needed, in 
which case one should obtain poly homogeneous expansions in which uniformity 
is lost when the corner r = x = is approached; this will be discussed elsewhere. 

Theorem 3.7 Consider a solution (<p,ip) G CqoXCqo of the system (3.5), suppose 
that in addition to (3.12), (3.13), (3.19), and (3.50a) we have 

B lu B-,E£,l,& G L°°{n xo , T ) , (3.53a) 



x=0 

,El-E 



id , = (d T - 2d x ) <g> id , (3.53b) 



x=0 



E%-E% 



x=0 ' ^ ± 



G x 1+<5 ^(^ , T ) , (3.53c) 

x=0 

E A G x*/£{n X0 , T ) • (3.53d) 



If 



a, b G x^ k 5 (n x0tT ) , v?(0) G x^ k 5 (M X0 ) , 

with G R, and if the initial data satisfy corner conditions to arbitrary order, in 
the sense that 

V i G N 0^(0), 4V>(0) G «£(M Xo ) , (3.54) 
for some (i-independent) A G R, then 

XQ,1 ) ■ 

In particular, if k = oo then the solution is polyhomogeneous. 

Remark: The class of initial data satisfying corner conditions to arbitrary 
order is rather large; for example, if an initial data set (<£>(0), ip(0)) satisfies 
them, and if /, g are arbitrary functions smooth up to boundary on the initial 
data hypersurface, then (ip(0) + f,ip(0) + g) will also satisfy those conditions. 
More generally, large classes of such initial data can be constructed using a 
polyhomogeneous generalization of the Borel summation lemma. 

Proof: The hypothesis (3.54) with i = and Proposition 3.1 show that for all 
r G [0, T] we have 

y>(r), V(r) G tf*(M xo/2 ) . 

Proposition 3.3 shows then that the hypotheses of Corollary 3.5 are satisfied, 
and the result follows. □ 
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4 The semi- linear scalar wave equation 

Let / be a solution of the semi-linear wave equation 

a s f = H(x^f), (4.1) 
here D fl is the d' Alembertian associated with g. Set 

(n-l) 

f = V- — f; (4.2) 
Letting g = S7 2 g as in (2.1), from (2.3) we obtain 

□fi/ = " §)f + n-^H(x^, . (4.3) 

Let g = r\ be the Minkowski metric; under the conformal transformation (2.4) 
one obtains from (2.5) that g is again the Minkowski metric, and (4.3) becomes 

□„/ = n-^H(x^, n^f) . (4.4) 

We shall assume that the initial data for / are given on a hypersurface £ C 
which, in a neighbourhood & of is given by the equation 

Sn^ = {y° = i}. (4.5) 

This corresponds to a hyperboloid in ^# given by the equation x° + 1 = 
Vl + a? • It is convenient to introduce the following coordinate system (x, v, r) 
in a ^-neighbourhood of J? + : 

r = y°-l/2>0, 

* = (E(y l ) 2 )"-y ^ ' 

n*(-u) G S' n_1 , with v = (z/ 4 ) denoting spherical coordinates on S" 1-1 . Equa- 
tion (2.5) gives 

n = x(2t + x + 1) « x . (4.7) 
If we let /i denote the unit round metric on S 1 " -1 , we then have 

7] = 2dxdT + dx 2 + (x + r + l/2) 2 h , (4.8) 

and 
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where is the Laplace-Beltrami operator of the metric h. We set 

e_ = d T , e+ = d T - 2d x , e A = {x+1 } +1 / 2) h A , (4.10) 

0_ = e_(/~), + = e+(/), (4.11) 

= Y>A = (a+r + 1/2) M/) , (4-12) 

where Ha denotes an /i-orthonormal frame on S n ~ l . We use the symbol D to 
denote the covariant derivative operator associated to the metric h. (The use- 
fulness of introducing two different objects for /i^(/)/(x + r + 1/2) will become 
clear shortly.) Equation (4.4) implies the following set of equations: 



(4.13) 
(4.14) 



e-fa+) -D e J) A - 2{X + T + 1/2 ) <I>+ = - 2 (x+r+i/2) ^- + Q + 
-e A ((f> + ) +e+(ipA) - ( a+T +i/2) ^ = & A 

e-(^) -eA (</>-) + ( g+T +i/2) ^ = oa , 

— D eA (j)A +e+(0_) + 2(^+1/2) ^- = 2(x+r+l/2) ^+ + 6 ~ ' 

e_(/) = </>_, (4.15) 

e+(f) = <!>+, (4.16) 

with oa = &a = and 

a+ = 6_ = -G= -r^ff^.nV/) . (4.17) 

4.1 Existence of solutions, space derivatives estimates 

We note that the partial differential operator standing on the left-hand-side of 
(4.13) is symmetric hyperbolic; the same holds true for (4.14), or for the joint 
system (4. 13)- (4. 16). Now, part of our technique consists in deriving weighted 
energy estimates for symmetric hyperbolic systems having the structure above, 
cf. Section 3. Each such system comes with his own estimates, so that for the 
systems (4.13) and (4.14) we can obtain estimates with different weights. This 
allows us to handle a reasonably wide range of non-linearities, giving existence 
and blow-up control for initial data in weighted Sobolev spaces (with conormal- 
type blow-up at 

Theorem 4.1 Consider Equation (4.1) on MJ 1 ' 1 with initial data given on a hy- 
perboloid 5? D S Xo ,0 in Minkowski space-time, and satisfying 

~ n — l 

/k , = O-— /|e xo ,o e ^(2x0,0), (4.18) 
d x (^^f)K,o e % a (^ xo ,o)n^ k a ~ 1/2 (^o), (4.19) 

n — l 

d T (n-— /)| E(t0i0 e ^ a (S XOi o), (4.20) 

with some fe > ^ + 1, — 1 < a < —1/2. Suppose further that H has a uniform 
zero of order I at / = 0, in the sense of (A. 30), with 

(-4 , n = 2 , 

£ > J 3 , n = 3 , (4.21) 

I 2 , n > 4 . 

Then: 
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1. There exists < r + < T (< xq/2), depending only upon x$ and a bound 
on the norms of the initial data in the spaces appearing in Equations (4.18)- 
(4.20), and a solution / of Equation (4.1), defined on a set containing Q xo ,t+ . 
satisfying the given initial conditions, and satisfying 

\\f\\L°°(Cl a0lT+ ) < 00 • 

2. Further, if r* is such that / exists on fi X0)Tt and satisfies ||/||x,oo(q t j < oo, 
then for < r < r* we have 

d T j\^ T G ^(£ X0 , T ) , dJ\^ T G J< _l (Exo,r) n ^(E XOiT ) , 
with a r-independent bound on all the norms. 
Remarks: 

1. Integration in x of condition (4.19) implies that / G L°°(£ XOi o)- 

2. Some further information can be found in Theorem 4.3 below. 

3. If the inequality in (4.21) is not an equality for n = 2, 3 (no further re- 
strictions for n > 4), then a proof similar, but simpler, basing on Proposition 3.2 
instead of 3.1, leads to the same result with a = —1/2. 

Proof: As before, we write ||/(r)||^a for ||/|e xo , t || j>t«(e xo>t ), etc. Recall that 
the standard theory of hyperbolic systems (cf., e.g., [35, Chapter 16, Vol. Ill] 5 ) 
shows that for any < x\ < xq there exists T(x\) > 0, satisfying 2x\ + T < 
xq, and a solution / of (4.4), defined on SI Xi>Xo ,t, with initial data on S XljXo 
obtained from those on £ xo by restriction. The idea of the proof is to derive x\- 
independent, weighted a priori estimates for the solution. These estimates will 
guarantee that the existence time T(x\) does not shrink to zero as x\ goes to 
zero; they will also guarantee that the weighted Sobolev regularity is preserved 
by evolution. We start with the following: 

Lemma 4.2 Under the hypotheses of Theorem 4.1, consider on fi xi xo t the- sys- 
tem (4.12)-(4.16), set 

E a (t) = \\fm 2 je? + U-(t)\\ 2 je? 

+\\Mt)\\ 2 «_i+Ell^(*)llk«- ^ 4 - 22 ) 

Then there exists a xi-independent constant C such that 

E a (t) < C ^E a (0)e ct + j\ c ^S(s)ds^ , (4.23) 

where 

S{s) = El|aA(s)||^ +||a + (s)||2 1/2 

A k 

+ ll&-Wll^a-l/ 2 +Ell^( a )H^- 1 • ^ 
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Proof: We wish, first, to apply Proposition 3.1 to the system consisting of 
Equation (4.14) together with e_(/) = <f>~; in order to do this we set 

*-(#')• *-*- 

We choose E^d^ = e± ® Id, we set 

(4 - 25) 

and we define 

E a {t) = \\f{t)fje k ° + \\e-(f)(t)\\ 2 ^ + \Mf)(t)\\*r k ° ■ 

A 

The hypotheses ^T — ^5 of Proposition 3.1 are readily verified, and for any 
a < — \ the inequality (3.20) gives 



I w . 
A 



£«(<) < C^E a (0)e ct + J^ e c ^ |M* 

+ U + (s)f Jifa _ 1/2 + HMs) ll;U-i/2) ds} • (4-26) 
Next, we apply Proposition 3.1 directly to (4.13): setting 

E a ,(t) = ||e+(/)(t)||^+X;il^(/)(t)||; 
for any a' < —1/2 it follows from (3.20) that 

E a '(t) < C^E a/ (0)e ct + J\ c(t '^{\\a + (s 

+ H-(s)\\% a ,- 1/2 + £ \\bA(s)\\% a ,- 1/2 ) cfa J . (4.27) 

We set 

E(t) = E a {t) + E a _ 1/2 {t) . 
It follows from (4.26) and (4.27) with a' = a - 1/2 that we have 

E(t) < C ^E{0)e ct + jf * e c (*- s ) (£(s) + 5(s))ds^ , (4.28) 

with S(s) as in (4.24). Equation (4.23) with E a replaced 8 by E follows now 
from GronwalPs Lemma. Since E a is equivalent to E, our claims follow. □ 

Returning to the proof of Theorem 4.1, Lemma 4.2 applied to (4.13)-(4.15) 
gives (recall that G was defined in (4.17)) 

E a (t) < C^E a (0)e ct + J*e c ^\\G(s)\\ 2 K _ 1/2 ds^ . (4.29) 



3 The constant C in Equation (4.23) does not necessarily coincide with that in (4.28). 



25 



By hypothesis the function H appearing in (4.1) has a uniform zero of order 
I > 2, in the sense of (A. 30); we wish to use (A.31^to control the term containing 
G(s) in (4.29). This requires an L°° bound on /, which will be obtained next. 
As k > n/2 + 1, the Sobolev embedding (A. 24) gives 



||e_(/)(s)|& f + |M/)( S )||; rV2 + \\e A (f)(s)\\^ a < CE a {s). (4.30) 
Now the conditions (4.21) on n and I give 

|G(r)| < C||/(r)||l. x^~ n -¥ < C\\f(r)\f Lac X ~W , 
so that (recall that a < — ^) 

||G(r)||^.<C||/(T)||lco. (4.31) 

From (4.13) we have 



d T <t> + -- — - — - — tt4>+ = D e .ipA 7 — - — — r -r<t>--G, (4-32) 

2{x + t + l/2y + 2(rr + r + l/2) y 



and (4.31) together with Proposition B.l yield 

U + {t)\w<> < Ce ct \\M0)\h o « 

+C f e C ^- s \\\D e Ms)\\^ + W-{s)\\ v « + \\G(s)\\ v «) ds 
Jo 

< Ce ct \\<P + (0)\\ K + f e c ^C{E a {s),\\j{s)\\ L ~)ds , (4.33) 
J o 

for some continuous function C(E a (-), Integration of d x f = \{4>- — 

cj) + ) over [x,xq — 2r] gives 

~ ~ 1 rxo—2r 

|/(r,x)|<|/(r,Xo-2r)| + -||(0_-^ + )(r)||^ J s a ds . 

For any < r < r* < xo/2 the /(r, xq — 2r) term is estimated by a multiple 
of the initial energy in a standard way, which leads to the estimate (recall that 
a > -1) 

||/(t)||l~ < CE a (T) + Ce CT U + (0)\\^ 

+ [ T e c ^C(E a ( S ),\\f(s)\\ L ~)ds . (4.34) 
J o 

Next, 



n-l 



IIGWII^a-i/a < C\\H(S,;X— f)\\^ a _ 1/2+ 



n+'3 , 



and our hypothesis that H has a uniform zero of order £ together with (A. 31) 
gives 



\\G(s)\\^ rl/2 <c(\\f(s)\\ L ~) 
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In view of (4.34) this can be estimated by a function of E a (s) and of ||/(s)||l°°, 



\\G(s)\\^ rl/2 < C(\\f(s)\\ L ~)\\f(s)\\^ ka 

< C (\\f(8)\\ L ~) E a (s) , (4.35) 

provided that 

I > ^ (4.36) 
n — 1 

(which coincides again with (4.21)). If (4.36) holds, from (4.29) and (4.34) we 
obtain 

||/(r)|U~ + E a (r) < Ce CT (s Q (0) + R/(0)||^ + R/(0)||^) 

+ jT $ (t, s, \\f(s)\\ L ~,E a ( 8 )) ds , (4.37) 

for some constant C, and for a function $ which is bounded on bounded sets. It 
then easily follows that there exists a time r+ and a constant M, depending only 
upon xq and a bound on the norms of the initial data in the spaces appearing 
in Equations (4.18)-(4.20), such that ||/(t)||z,°° and E a {r) remain bounded by 
M for < r < r+. Since all the objects above were xi-independent, so is r + . 
By the usual continuation criterion (c/., e.g., [35, Proposition 1.5, Chapter 16, 
Vol. Ill] 9 ) the solution exists on fi XljXOjT+ for all x\; it thus follows that the 
maximally extended solution of the initial value problem considered here exists 
on a set which includes Q Xo ^ T+ . 

To establish point 2, suppose that a global a-priori L°° bound on / is known. 
Then (4.29) and (4.35) give a linear integral inequality on E a , and GronwalPs 
Lemma gives a global bound on E a . Arguments of the last part of the proof of 
point 1 yield the result. □ 

For the purpose of estimating time derivatives of the solutions we will need 
a generalization of Theorem 4.1, which covers the equations contained by time- 
differentiating Equations (4.13)-(4.16). There are lots of ways to relax those 
hypotheses; for simplicity we shall only make those generalizations which are 
strictly necessary for the arguments in the next section to go through. First, the 
fact that / is scalar valued plays no role in our considerations above; henceforth 
we assume that / has values in for some N > 1. Next, the definitions (4.10) 
of e± and e\ will be kept. We will consider systems of the form 





B12 




B21 


B22 


)(; 


( 0+ 

V 4>A 


)• 





f=[ Y. ) , tf={ J" ) (4.38b) 



9 In that reference symmetric hyperbolic systems on a torus are considered; however simple 
domain of dependence considerations show that the results there apply to the setup here. 



27 



together with 



^ = (3!+T +i/2) M/) + b aJ , (4.39a) 

^ = {x+T l +1/2) h A (f) + B A J , (4.39b) 

e_(/) =B cf>_ + Bj, (4.39c) 

e+(7) = 0+, (4.39d) 



for some matrix valued functions i?^, £?o, -Bi, with Bq — invertible. 
Here 

e_ £ A D A 

(£ A yD A e+ 



(4.40) 



is the (geometric) principal part of Equations (4.13)-(4.14). The nonlinear term 
G = G(x", f) will be labeled as 

G = (G ! e+ (^_),G e+ (^ A ),G e _(^ A ),G e+ (^_)) , (4-41) 

with the order of the components following that of Equations (4.13)-(4.14). The 
Bab's will be labeled as B ( f > _ t( f >+ , B ( f ) ^ A , etc.; for example, in this notation, the 
second of Equations (4.14) takes the form 

e + (<f>-) = D eA 4>A ~ B (j> _^_(j)- 

-#4>_,4> + 0+ - B^_^J> A - B 4> _^ A '4) A + 6_ + G e+ (0_) , (4.42) 

with actually = B^_^ A = 0. 

Some effort will be needed to prove the information of point 3 of the theorem 
that follows; this is needed to be able to iteratively apply that theorem in the 
next section: 

Theorem 4.3 Consider the system (4.38)-(4.39) with 

IKt)H^ + \\b(T)\\^ + sup \\B ab (T)\y>+ sup \\B A: x{r)yo 

a,6=l,2 k A=l,2; \=<j>,ip k 

+ ||B (r)||^o + \\B q \t)\\ l ~ + WB^t)^ < C , (4.43) 

for some constant C, and suppose that 

G(x», /) = n-( n+z )/ 2 H (x», n^-V/ 2 /) , (4.44) 

with G e _(^ A ) = 0, and with H having a uniform zero of order £ in the sense of 
(A. 30), with £ satisfying (4.21). If the initial data satisfy (4.18)-(4.20) with some 

fc>§ + l,-l<a:< -1/2, then: 

1. The conclusions of point 1. of Theorem 4.1 hold with a time r + depending 
only upon the constant C in (4.43) and a bound on the norms of the initial 
data in the spaces appearing in Equations (4.18)-(4.20). 

2. The conclusions of point 2. of Theorem 4.1 hold. 
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3. Under the hypotheses of point 2. of Theorem 4.1 we also have 

||(a: + 2r)a T 7|| L =o (nisoiT . ) <oo. (4.45) 

Remarks: 1. The condition G e _(<^ A ) = can be weakened to 

G e _ iM (x»J) = n-^ +2 )l 2 H e _ {4>A) {x^M n - 1),2 J) , (4.46) 

for some function H e _^ A ^ with a uniform zero of order i. Similarly it suffices 
to assume that 

G e+ ^ A) {x»J) = n-^y 2 H e+{M (x»,rt n - 1)/2 f) , (4.47) 

for some function H e + (if, A ) with a uniform zero of order L 

2. If the inequality in (4.21) is not an equality for n = 2, 3 (no further 
restrictions for n > 4), then the result remains true with a = — 1/2, see Remark 
3. after Theorem 4.1. 

Proof: Let us start by remarking that, because tpA = 4>A, in equations such 
as (4.42) we can replace B < j > _^ A by B ( j > _ t< j >A + A obtaining a system in 

which B ( f ) _^ A = 0. Proceeding similarly with the other equations we may thus 
without loss of generality assume that 

B^ A = . (4.48) 

The proof of points 1 and 2 is then identical to that of Theorem 4.1, with the 
following minor changes: Equation (4.32) is replaced by the equation 

e- (</>+) +#</>+,</>+</>+ = 

D eA <j) A - B^_^_<\>- - B^ + ^ A 4> A + a+ + G e _ {(j)+) (4.49) 

to which Proposition B.l still applies, recovering (4.33). Further, the equation 
d x f = (</>- — </>+)/2 has to be replaced by 

~ Bi~ B (p_ - <f> + 
dxf + Y f = 2 ' 

and the desired conclusion is obtained by Proposition B.3. The remaining 
arguments do not require any modifications. 
To prove point 3, from (4.42) we obtain 

e+[(x + 2r)<f>-] = 

(x + 2r) (D eA 4> A ~ B^,_^_<j)- 

-B^_^J>a~ B 4> _^ + 4> + + b- + G e+ (0_)) , (4.50) 
From Equations (4.31), (4.33), and (4.39c) together with 

<f>-,(j> A € W£ C ^ , D eA 4> A € J^L, C ^ , 

we obtain 

e + [(x + 2r)^_] < Cx~ a , 
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for some constant C depending only upon the initial data and ||/||i,oo(n t j. 
Integrating as in the identity (3.41) we arrive at 



IB; 



3 1 {(x + 2T)(d T f-Bj)(x,v,T)}\ 

< {B, 1 {(x + 2r)dJ(x + 2t,v,0)} \+C (||/|| L -(^ , T J + c) 

< c (\\d T f\\^ + \\fm L ^ XOiT j + c) , 

and Equation (4.45) follows. □ 



4.2 Estimates on the time derivatives of the solutions 

So far we have established existence of solutions with initial data in weighted 
Sobolev spaces, as well as weighted estimates on the space-derivatives of the 
solutions. The next step in proving polyhomogeneity is to establish estimates on 
time-derivatives. Similarly to the linear case, the question of corner conditions 
arises. In order to handle that, we introduce an index m, which corresponds to 
the number — perhaps zero — of corner conditions which are satisfied by the 
initial data. Next, the definition (A. 30) of a uniform zero of order I has to be 
strengthened by adding conditions on time-derivatives: we shall require that 
for all MGl, \p\ < M, < i < min(fc, I) and for all < j < m there exists a 
constant C = C(M, m, k) such that we have 



d^F(r,;p) 



dp l dri 

We start with the following: 



< C\p\ 1 - 1 . (4.51) 

f!L ■ ■ 



Theorem 4.4 Let N B m > 0, consider a solution / : fix ,r» ^ K of Equation 
(4.1) satisfying 

Wfh°°(n X0 , T ,) < oo , 

and suppose that 

0<i<m + l 4/| Exoj0 g ^Wi-i^o) , (4.52) 

o<i<m a^/ls^o e ^(s IO , )n4 Q + -^(s XO)0 ), (4.53) 

with some k > ^ + 1 and — 1 < a < —1/2. Suppose, further, that H is smooth in 
/ and has a uniform zero of order t at / = 0, in the sense of (4.51), with i as in 
Equation (4.21). Then for < r < r* and for < i < m, < j + i < k + m — n/2 
we have 

[(r + 2x)d T ydif\^ T € L°°(E XOiT ) n J^m+i-i-ji^r) , (4.54a) 
d x [(r + 2x)d T Mf\v a0tT G jCli-i( E *o,r) n ^(E, , r ) , (4.54b) 

and 

< p < k - n/2 [(r + 2x)a T ]^ T m+1 7|E :co , T G ^_ P (S, ,r) , (4.55) 
with r-independent bounds on the norms. 
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Remark: As before, in dimensions n > 4 the result remains valid for a = —1/2; 
in dimensions n = 2, 3 the value —1/2 for a is allowed if the inequality in (4.21) 
is not an equality. 

The proof below actually proves the analogous result for the systems con- 
sidered in Theorem 4.3, provided that obvious time-derivative conditions on the 
coefficients are added to (4.43), the simplest possibility being 

a, 0=1, 2 

+ sup ||4S A , A (r)|U . + ||4Bo(t)|Uo . + ||4Bi(t)|Uo . < 

(4 

with < % < m + k; the same remark applies to Corollary 4.5 below. Before 
passing to that proof, we note that an important consequence of Theorem 4.4 
is that corner conditions will hold at any time r > 0, regardless of whether or 
not they hold at r = 0: 

Corollary 4.5 Under the conditions of point 2 of Theorem 4.1, for any < r < r* 

and for < z < fc — 1 — n/2 we have 

4/k , T g L°°(Z X0 , T ) n ^ 1 _ i (£ X0)T ) , 
a^/k , r G <7(E IOlT ) n«3*(E XOiT ) . 

We shall need the following simple Lemma: 

Lemma 4.6 Let F(x fl ,p) be a function which is smooth in p at fixed x^ and 
suppose that it has a uniform zero of order I > 1 in p. Then 

1. For all i G N the function d\ (F(x ,x , u(x^)) has a uniform zero of order £, 
when viewed as a function of (u, d T u, . . . , d\u). 

2. Let # = d p F, then has a uniform zero of order £ — 1. 

Proof: Let u = smoothness of -F in p allows us to write 

F(x,T,u) = A il ... it u il ...u it , (4.57) 

with some coefficients -4^...^ = -Aji...^ t, u) which are smooth in u, and totally 
symmetric in ii, . . . recall that the summation convention is used through- 
out. Point 2 immediately follows from (4.57). From that equation we also 
obtain 

d T F(r, x, u) = (d T A il ... ii + d u iA il ... it d T u l ) u n . . . u H 
+£A h ... k u h ...u k -^d T u k , 

which proves point 1 for i = 1. The result then follows by straightforward 
induction. □ 
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We can pass now to the proof of Theorem 4.4: 

Proof: We assume that Equations (4.38)-(4.39) are satisfied; Theorem 4.3 
shows that (4.54)-(4.55) hold with i = j = p = 0. Consider the vector-valued 
function 

(/, {x + 2r)d T f, <p, {x + 2 T )d T ip, Y>, (x + 2r)a T ^) , 

so that the new function / in (4.39) is (/, (x + 2r)d T f), while the new functions 
ip, respectively ij>, in (4.38b) are (if, (x + 2r)d T (p), respectively (ip, (x + 2r)d T ip). 
We claim that a set of equations of the form (4.38)-(4.39) holds for those new 
functions. Consider, for instance, Equation (4.39c); set 

f:=(x + 2T)d T f, 0- := (x + 2r)d T <\>- , 

etc., we have 

e-(/) = d T ([x + 2T){B ^ + Bj)) 

= B o 0- + (2B + (x + 2r)a TJ B o )0- 
+BJ+ (2B 1 + (x + 2r)a r B 1 )/, 

which is linear in (/, /, </>_, 0_). In fact 



e_ 




+ 



B 

2B + (x + 2T)d T B B 

B x 

2B X + (x + 2T)d T B 1 Bi 




and the new matrix Bq is again invertible, as desired. Next, 

e_(0+) = d T {{x + 2T)d T <P + ) 

= d T ((x + 2T)(-D A <j> A -B 4>+( p_<f>_ 

-B<l> + <p A 4>A - #</>+</>+</>+ + a + + G e _(0 + ))) 
= —Da4>a — B < f >+( f ) _4>^ — B^^^a — B < f >+( f )+ (j) + 
+ linear (99,^) + a+ + G e _$ + ) , 
a+ = -2L» A 0A + <9 T a+ G ^fc+ m _i , 
G e_(^ + ) = a T (G e _ (0+) )(x + 2r), 

where "linear" denotes terms which are linear in the relevant variables. The 
equation for e_ ((Pa) is handled in a similar way. The equations involving only e+ 
or &a are straightforward, since those operators commute with multiplication by 
(x + 2t) . By Lemma 4.6 the new non-linearity has again a zero of order £, when 
considered as a function of (/, (x + 2r)d T f). In order to apply Theorem 4.3 we 
need to check that the initial data are in the right spaces. Clearly 



((* + 2r)0 T /)(o) = xd T f(o) g jt^ 1 c je k \ m n l 



(d x {{x + 2r)d T f)) (0) = (d T f + xd x d T f) (0) G ^ +m _! C n Jt£ 



a-1/2 
m—1 
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Condition (4.20) requires some more work: 



(d T ((x + 2r)d T fj)(0) = (2dJ+ xd 2 j) (0) 

= (2dJ+ x{2d x + e + )d T f) (0) 

= (2dJ+ 2xd x d T f + xe + (b Q (/>- + (0) . 

The first two terms are obviously in J^ i _ m _ 1 , and so is xe + (B\f) = x(<9 T — 
2d x )(B-f). Equation (4.42) gives 

(xe + (4>- )) (0) = x (£> eA A - -B<^_ ,</>_ (t>- - B 4> _ )4>+ <f> + 

-B 4> _^ A <j)A - B^^t/ja + b- + G e+{0 _ ) ) (0) . 

The desired property (xe+ (£>()</>_)) (0) € ^fc+ m _i follows immediately; the only 
non-trivial term is xG e+ ^_y the ^,+ m+1 norm of which can be estimated by a 
function of 1 1 /(O) 1 1 L oo and ||/(0)||.*>« , cf. Equation (4.35). Now, (x + 2r)d T / 
is uniformly bounded on £1 Xo ,t, by point 3 of Theorem 4.3, so that we can apply 
point 2 of Theorem 4.3 to conclude that Equations (4.54)-(4.55) hold with 
j = p = 1 and m = 0; straightforward induction establishes Theorem 4.4 for 
the remaining j's and p's. 

Consider, now, m = 1; the result already established with m = shows 
that d T f(r) exists and satisfies (4.54) with i = 1 for any r > 0; similarly (4.55) 
holds with m = 1 for any r > 0. JNTow, a calculation similar (but simpler) to 
the one done above shows that (f,d T f) satisfies a system of equations of the 
form (4.38)-(4.39) with initial data satisfying the conditions of Theorem 4.3 by 
hypothesis; the uniform bounds on some interval [0, r+) follow by point 1 of 
that theorem. We therefore have 

||(/,3 T /)||Loo (n<S0iT ,) < 00 • 

We can then apply the result already established for m = to the system of 
equations satisfied by (/, d T f) to obtain the conclusion of Theorem 4.4 with 
m = 1. An induction upon m finishes the proof. □ 



4.3 Polyhomogeneous solutions 

The aim of this section is to establish poly homogeneity of solutions of a large 
class of semi-linear systems of the form 

d T (p + Bmp + B 12 ip = L U ip + Li 2 Tp + a + G ip , (4.58a) 
d x tp + B 2 i<p + B 22 ip = L 21 (p + L 22 tp + b + G^ , (4.58b) 

with a nonlinearity 

G = (Gip, G^) 

of the form 

G = x^Hix", x q ^ l ,x qS+1 ij 2 , x" 5+ V) . (4.59) 
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Here we have decomposed tp as 

*-{*,)■• (460) 

this is motivated by different a priori estimates we have at our disposal for the 
appropriately defined components ipi and ip2 of ip in the applications we have 
in mind. Polyhomogeneity of solutions of (4.1) will follow as a special case, 
see Theorem 4.10 below. We will need to impose various restrictions on the 
function H, in order to do that some terminology will be needed. We shall 
say that a function H(x^,u) is 5-polyhomogeneous in x with a uniform zero of 
order I in u if H is smooth in u G at fixed x^, if H satisfies (A. 30) for any 
< i < min{/, k} and any k G N, if 

ViGN diH(-,u)e^i 

at fixed constant u, and if we have the uniform estimate for constant u's 

Ve > 0,M > 0,i,k G N 3 C(e, M,i,k) V|u| < M < C(e,M,i,k) . 

(4.62) 

The qualification "in n" in "uniform zero of order Z in n" will often be omitted. 
The small parameter e has been introduced above to take into account the 
possible logarithmic blow-up of functions in g/^ at x = 0; for the applications 
to the nonlinear scalar wave equation or to the wave map equation on Minkowski 
space-time, the alternative simpler requirement would actually suffice: 

VM>0,i,^H 3 C(M, i, k) \/\u\ < M \\d^H(-,u)\\^o < C(M,i,k) , 

(4-63) 

again for constant u's. Clearly functions which are jointly smooth in u and 
in x^ satisfy the above conditions; Lemma 4.7 below provides another class of 
such functions. The following simple facts about functions in the above class 
will be useful: 

Lemma 4.7 Let mi,m2,k G N, m\ < rri2, and let P(x^,u) be a polynomial in 
u = (u 1 ,... u N ) of the form 

P(x",u)= Ph...i j {x lx )u il ...ifi , 

mi <j<rri2 

with coefficients Pi 1 ...i j {x^ 1 ) G si^. Then: 

1. P is (5-polyhomogeneous in x with a uniform zero of order mi. 

2. If 

for some A > 0, then for any e > we have 

P(.,x^/)ex m "V^<" £ ). 

The proof of Lemma 4.7 is elementary and will be left to the reader. 



(4.61) 
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Lemma 4.8 Let k,q G N and let H(x fl ,u) be 5-polyhomogeneous with respect 
to x with a zero of order m in u. If 



n l°° + , g = o, 

otherwise, 



for some A > 0, then for any e > 

ff(-,//)Gx^(^ + ^- £ ). 
Proof: We Taylor-expand in u to order r, where r is any number satisfying 

rq8 > mq8 + A . 

We then have 

H(x»,x qS f) = P(x»,x qS f) + R , 

where P is a polynomial and R is a remainder. We note that the coefficients 
of the expansion of P can be obtained by differentiating with respect to u and 
setting u = 0, and are therefore in =e^, by (4.61). Further, the usual integral 
formula for the remainder in a Taylor expansion together with (4.62) shows that 
R has a uniform zero of order r, in the sense of Equation (A. 30). The result 
follows from Lemma 4.7 and from Lemma A. 5. □ 

We are ready now to pass to the proof of the non-linear analogue of Theo- 
rem 3.4: 

Theorem 4.9 Let p G Z, q, 1/5 G N, -1< (3 f G R, k G N U {oo}, and let 

e tf£(nx ,T) x v£(n XOtT ) , ^ g l°°(si xo>t ) 

(ipi as in Equation (4.60)), be a solution of (4.58) with G of the form (4.59), where 
H is (5-polyhomogeneous in x with a uniform zero of order 



r> — I 

m > £ . (4.64) 



9 

Suppose that Equations (3.43)-(3.44) hold, and that 

B u G « n L°°) (fi X0)T ) , B 12 , 5 22 , B 21 G <(^ ,t) , (4.65a) 
a, 6 G <(^ ,t) , ¥>(0) G <(M Xo ) . (4.65b) 

Then 

V G (x^-^V/ + ^f) (^ , T ) = x min((mp-,)«,0)^ ) 

^ G x mi„ {(m ^p) 5+ i,i}^ ( ^ oiT) + Coo (TQ c n (n X0)T ) . 

If one further assumes 

L5 t 2 ,B 12 ,a,^(0),G^(-,0)GL oo (^ , T ), 
then it also holds that 

<p G ( x (^-p) V/ + n (n X0)T ) . 
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Remark: Obviously the theorem remains true if we replace G by a finite sum 
of nonlinearities satisfying the above hypotheses, with different p's and g's for 
each term of the sum. 

Proof: The result is established by a repetition of the proof of Theorem 3.4, 
using Lemma A. 5 and Lemma 4.8 to obtain the necessary estimates on the 
non-linear terms. We simply note that the condition on the order m of the 
non-linearity guarantees, using Lemma A. 5, that 

d x j> = c 2 e ^ , 

with 

A = mmf/?', mq8 — p5} > —1 , 

hence ip G L°° by integration. Decreasing (3' if necessary we may without loss 
of generality assume that /?' = A. When applying Lemma 4.8 it is convenient 
to view the function H as a function of the variable / := (ipi,xip2,x(p) G L°°. 
The remaining details are left to the reader. □ 

As a straightforward corollary of Theorem 4.9 one obtains: 

Theorem 4.10 Let S = 1 in odd space dimensions, and let 5 = 1/2 in even space 
dimensions. Consider Equation (4.1) on R"' 1 , n > 2, with initial data 

7| {r=0 } , df/dr\ {T=0} G n (M xo ) . 

Suppose further that H{x^, /) is smooth in / at fixed x M , bounded and 5-polyhomogeneous 
in at constant /, and has a zero of order £ at / = 0, with i as in (4.21). Then: 

1. There exists r + > such that / exists il XOtT+ , with 

ll/IU»(n. ,r + ) ■ ( 4 - 66 ) 

2. If the initial data are compatible polyhomogeneous in the sense that there 
exists A < 1 such that 

Vi G N 3^/(0) G tf- A (M X0 ) , 

then the solution is polyhomogeneous on each neighbourhood Q xo ,t* of ^ + 
on which / exists and satisfies (4.66) with r + replaced by r*. 

Proof: Point 1 is Theorem 4.1 specialized to polyhomogeneous initial data. 
To prove point 2 we set 




(4.67) 



and 

if = 4>+ . (4.68) 



36 



Then Equation (4.3) takes the form (4.58) with 

G = -n-^H^, n^f) = -VL~^H{x^ : O^Vi) , (4-69) 
G v = -G, G^=0, G f2 = ( " Q G \ . (4.70) 

For n even we take 5 = 1/2, p = n + 3, q = n — 1; the condition (4.64) then 
reads m > which coincides with (4.21). For n odd we take 5 = 1, p = ^y^, 
q = and (4.21) guarantees again that (4.64) holds. □ 



5 Wave maps 

Let (^V,h) be a smooth Ricmannian manifold, and let / : (^#,g) — ► (jy,h) 
solve the wave map equation. We will be interested in maps / which have 
the property that / approaches a constant map /o as r tends to infinity along 
lightlike directions, fo(x) = Po S for all x G . Introducing normal 
coordinates around po we can write / = (f a ), a = 1, . . . , N = dim^, with the 
functions f a satisfying the set of equations 

Qfb Qfc 

□ s r + B^ c (/)^ I ^ 7 = 0, (5.1) 

where the T^'s are the Christoffel symbols of the metric h. Setting as before 
j a = n-^r, q = n 2 Q, we then have from (2.3), 



(5.2) 

In particular if (^#, g) is the Minkowski space-time (and if we use the same 
conformal transformation as in Section 2) we obtain a system of Equations 
(4.13)^(4.17) with = 6,4 = 0, with the obvious replacements associated with 
/ — > / a , and with G in (4.17) replaced by 

G a := -r^ 2 ?/) {^(-/ + «f_ + <^) 

- (n - l)n^f c \[x<j) b + - (1 + x + 2r)<^) - (n - l)/ 6 ] } . (5.3) 

5.1 Existence of solutions, space derivatives estimates 

As before, for even space-dimensions n the occurrence of non-integer powers of 
f2 above does not allow the use of the standard conformal method except for 
special target manifolds (^V,h), cf. [11]. This can be handled in our approach, 
and we show: 
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Theorem 5.1 Consider Equation (5.1) on W 1 ' 1 with initial data given on a hy- 
perboloid 5? D ^ xo ,o in Minkowski space-time, and satisfying 

r a| =n - a r 1 fa, e {feni»)(s xo ,„), n>3, 

/K,o-» / l^.o e \ (je« , n «f ) (E^o) , n = 2, (5 ' 4) 

77,-1 

9x(«— 5-/°)|e bo , e ^(S XOi o), (5.5) 

a ')l*Wo l(4«nL»)(S %0 ), n = 2. (5 - 6) 



for some fc > § + 1, -1 < a < -1/2. Then: 

1. There exists r + > and a solution / a of Equation (5.1), defined on a set 
containing fi X0)T+ , satisfying the given initial conditions, such that 

nrik° ( nw < 00 . n = 2 ( 5 - 7a ) 

r 

||are + (/ o )|| L oc (nit0iT+ ) + \\ xX if a h°°(n X0 ,r + ) 

i=l 

+ \\f a \\ L ~(n X0 , T+ ) + \\xd T f a \\ L ~ {Qxo!T+) < oc , n > 3 . (5.7b) 
Here the Xi's are the vector fields defined in Section A, cf. Equation (A. 7). 

2. Further, if r* is such that f a exists on £1 Xo ,t* with (5-7) holding with r + = r*, 
then for all < r < r* we have 

uniformly in r. If n = 2 we also have uniform bounds in the following spaces 

7 a | s , , T g « n jr fc %) (s X0)T ) , dJ a \^ T € n l°°) (e X0)T ) . 

Remark: Integration of condition (5.5) implies of course that / G L°°(E XOj o). 

Proof: The proof is similar to that of Theorem 4.1, but simpler, because we 
do not need to gain a 1/2 in the decay rate, as done in Lemma 4.2. We write 
Equation (5.1) in the form (4.12)-(4.16), with a a = &A = and with G in (4.17) 
replaced by G a defined in (5.3). We write G a as 

G a = A a + B a + C a + D a + E a , (5.8) 

with the order of terms in (5.8) corresponding to that in (5.3). Since we are 
working in normal coordinates, T? has a uniform zero of order one in the sense 
of (A. 30) at f a = 0. We want to use Equation (3.20) to get an a-priori estimate 
for the solutions of (5.1); for this we shall need to estimate the Jtfj^ norms of 
all the terms which occur in (5.8). The simplest such term is E a : 



(n - iy\\r a bc (n—f)(n—f)(n—f 



7b \ 



ot+(n+l)/2 , 



k 
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where we have used the fact that 0,/x is a smooth, and therefore bounded, func- 
tion. The function T^Q— /)(fi"r f c ){9.— f b ) can be viewed as a smooth 

n— 1 ' — ' 

function F of x~ z~ f a with a uniform zero of order three. We can thus apply 
(A. 31) with I = 3 to obtain 

||^)||^ < C(||/( S )|Uoc)||7||^a +2 -„ 

< C(||7( S )||^)||/||^« , (5.9) 

since n > 2. We note that in dimensions larger than or equal to three we 
have at least one power of x "left unused" above, which will be made use 
of in estimating the remaining contributions to G a . We proceed in a similar 
way with the other terms; in space dimension n = 2 we view 

n (n+l)/2 D a = 

n(n+i)/2( n _i^ 1 + x + 2r)n 21 T i T% c (n 11 ^ 1 f)f c 4> b _ as a smooth function F with 
a uniform zero of order three of {x^~ f a , x^~0" ), which leads to the estimate 

\\D(8)\\#« < C(||/( S )|| L oc,||0_( s )||^)(||/||^ + ||0_( s )||^) (5.10) 

On the other hand, in dimension 3 or higher we can view 

n (n+l)/2 D a 

as a 

function F with a uniform zero of order three of (x^~ f a , x~^~ x$ a _ ), which 
implies 

\\B{a)\\*~ < C(||/(a)||L<», ll^-(a)IUoc) (ll/H^a + ||x0_( S )||^)(5.11) 
Regardless of dimension we view f2( n+1 )/ 2 C a = ft( n+1 )/ 2 (n-l)xf^r£ c (ft^7)7 c <^ 

n— 1 ' ' n — 1 

as a smooth function with a uniform zero of order three of (x^~ f a , x^~ X(j)% ), 
obtaining thus 

< C(\\f(s)\\ L ~, ||^+( a )|| L =o) (\\f\U k « + ||z<M*)ll^r)(5.12) 

Viewing B a as a function of (x^2~ f a , x^~xcf>'\), and viewing A a as a function 
of (i~f, x~5~x(/) a L , x^~X(/)" ), one similarly obtains for n > 3 

POOH^* < C(||7(s)||loc,||^_( s )|| L oo,||^ + ( s )|| L oc)x 

(ll7ll^ + ll^-WII^- + lk<Ms)lk«) , (5-13) 

< C(\\f{s)U~, \\x<f>A(s)h~) (WfWsee + ||^( S )||^«)(5.14) 
while in dimension 2 it holds that 

\\A(s)\U« < C7(||7(s)||l<»,||^_(5)|Uoc,||^ + ( S )|| l =c)x 

(\\f\U? + U-( a )\U k a + • (5.15) 

\\B( S )\\^ < C(||7(a)||L=oJ|^(5)|Uoc)(||7||^« + ||^(a)||^) (5.16) 
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Summarizing, in space dimension two we have obtained 

\\G(8)\Ua < C(||/( a )|| L =o, \\</>-(a)\\ LOO , Ua(s)\\l°°, ||x^ + ( a )|| L oc) x 
(ll/ll^* + \\<P-( s )\\j? k a + \\x<t>+(s)\\jf? k a + Ua{s)\\j%>«} 
< C(||/( a )|| L =o, U-(8)\\ Loo , Ua(s)\\lo°, x 

V^a(s) , (5.17) 

where 

E a it) = mm^ + H-m 2 ^ 

+\\Mm 2 ^+J2\\^\\h a ■ (5 - 18) 

A 

On the other hand in higher dimensions we can write 

||G(s)||^a < C(\\f(s)\\ L oo,\\x(t)A{s)\\ L oo,\\x(j)-(s)\\ L oo,\\x(l) + (s)\\ L oo) x 

y/E a (s) . (5.19) 

To obtain a closed inequality from Equations (3.20) and (5.17) or (5.19), we 
need to control all the norms occurring there. Since k > n/2 + 1, from 
Equation (5.17) and the weighted Sobolev embeddings we obtain 

\\G(s)\\v° < C(||7( a )|| L =o,||^_( S )|| L oo, U A (s)\\ L oo,E a (s)), (5.20) 

in n = 2, or — from (5.19) — 

\\G{s)\U« < C(\\f(s)\\ L oo,E a (s)), (5.21) 

for n > 3. The identity 

f a (r,x) = f a (r,x - 2r) - - J {<f> a _ - <f> a + )(r, a ) ds (5.22) 

yields 

< c(V^(O) + ||0-^)lk o « + ||0 + (s)|k«) 

< C (yEjO) + y/E a (s)) . (5.23) 
for n > 3, while if n = 2 we use the estimate 

< C (V^a(0) + V^OO) • (5-24) 

In Equations (5.23)-(5.24), for notational simplicity we have estimated / a (r, xo — 
2r) and its angular derivatives by a multiple of the initial energy E a (0) ; strictly 
speaking, this should be some functional of (E a (0),T*) for r* small enough; then 
such an estimate holds by standard methods for < r < r* < xq/2. Further, 
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such an inequality is correct if we already have a weighted L°° bound as as- 
sumed in point 2. of the theorem. If n > 3 Equations (3.20) for a < —1/2 or 
(3.32) if a = -1/2, (5.21) and (5.23) give 

E a (r) < CE a (0) + f T <& (E a (s)) ds , (5.25) 
Jo 

for some constant C, and for a function <I> which is bounded on bounded sets, 
and we conclude as in the proof of Theorem 4.1. 
If n = 2, we note the identity 

(j)_(T,x) = </>_(0,x + 2r) + f e+ (<?!>_ )(ct, 2 (r-cr) + x) da . (5.26) 

Jo 

From the second of Equations (4.14) we obtain 

\e + (<t>_)(s,x)\ < C (||<M*)lk f + Ha(s)\\^ + ||<M*)lk" + \\G(s)\\^) x a , 
so that 

\<P-(t,x)\ < \\^(0)\\ LX+C £ (U-(a)y (? + \\Mv)\W? + \\<f>+(<7)y " 

+ \\G{a)\\^) {2{r - a) + x) a da. (5.27) 

It follows that 

HMt)||l=c < U-{m^+C (^fEja) + \\G{a)U«) (r - a) a #.28) 
Let 

= + IIM*)IU« + + V 7 ^) • (5.29) 

It follows from (3.20), (5.24) and (5.28) that we have 

F{t) < CF(0) + [ T $(F(a)) (1 + (r - <j) q ) da , (5.30) 
Jo 

where $ is a function bounded on bounded sets. We have the following: 

Lemma 5.2 There exists a time r*, depending only upon C , F(0), and the func- 
tion <I>, such that any positive continuous function F : [0, r+) — > M satisfying 
the inequality (5.30) with a > — 1 is bounded from above by C-F(O) + 1 on 
[0,max(T + ,T*)). 

Proof: Let 

M = sup \$(x)\ ; 

0<s<C*F(0)+l 

if M = the result is obviously true, so assume that M / 0. From Equa- 
tion (5.30) we obtain that on any interval [0, r) on which F < CF(0) + 1 we 
have 



F(t) < CF(0) + J M (1 + (r - ct) q ) da = CF(0) + M (t + 



T a+l s 

a + 1 
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(Equation (5.30) with r = shows that CF(0) > F(0), and continuity of F 
implies that the set of such intervals is non-empty.) The result is established 
by choosing 



= min 





"a + r 


l/(a+l)\ 


I 2M ' 


. 2M _ 





Because the existence time r* in Theorem 5.1 does not depend upon x\, 
Theorem 5.1 with n = 2 follows again by an argument identical to the one 
given at the end of Theorem 4.1. □ 

As in the case of the nonlinear wave equation (4.1), in order to obtain time 
derivative estimates we shall need a more general version of Theorem 5.1. Thus, 
we consider systems of the form (4.38)-(4.40) with a rather more general form 
of the non-linearity G appearing there. It should be clear from the proof of 
Theorem 5.1 that it is convenient to treat the case n = 2 separately, this will 
be considered in Section 5.3 below. We thus start with a result which holds in 
dimensions n > 3; the same proof gives similar results in dimension n = 2 for 
equations with a nonlinearity of higher order: 

Theorem 5.3 Let n > 3 and consider the system (4.38)-(4.39) with 

IK 1 ") 11^ + \\b(r)\\^ + sup ||B o6 (r)||^o 

a, 6=1, 2 

+ ||S (r)||^o + \\B^(t)\\ l ~ + ||Si(r)||^o < C , (5.31) 
for some constant C, with the nonlinearity G in Equation (4.38a) of the form 

(5.32) 

with G e _(<^ A ) = (cf. Equation (4.41)), and with H having a uniform zero of order 
t > 3 in the sense of (A. 30). Suppose that the initial data satisfy 

/1e, , ^ fr^r |e, , e feni°°)(M, (5.33) 
aj a |s, , e ^f(£*o,o), (5.34) 

d T f a k X0 , € Jf k a (X X0 , ) , (5.35) 

with some k > \ + 1, -1 < a < -1/2, then: 

1. There exists r + > 0, depending only upon the constant C in (5.31) and 
a bound on the norms of the initial data in the spaces appearing in Equa- 
tions (5.33)-(5.35), and a solution f a of Equations (4.38)-(4.39), defined on 
a set containing £1 Xo ,t + , satisfying the given initial conditions, such that 

r 

\\xe + {f a )\\ Lx{nxij T+) + \\xXj a \\ L ~ inxo T+) 

i=i 

+ \\f a h^n X0 , T+ ) + \\xd T r\\ L ^^ T+) < oo . (5.36) 
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2. Further, if r* is such that f a exists on ^I xo ,t, with (5.36) holding with r + = r*, 
then for all < r < r* we have 

/ a k ,. € L°°(E X0)T ) n J? k a +1 (X XOiT ) , (5.37a) 
Sr/ a |s, , T G ^(Sx 0lT ) , (5.37b) 
^7 a |s, , r G ^ Q (Sx 0) r) , (5.37c) 
with uniform bounds in r; this implies 

11x8^+ 1| L oo (nit0 iT .) + Hx^^IIloc^ t j + ||(x + 2T)d T f a \\ Lx(nxo ^ ) < oo . 

(5.38) 

If k > n/2 + 2 then we also have 

1 1 a; (a; + 2r)d T 4)- \\L°°(n X0 , Tt .) < 00 • (5-39) 

Proof: The transition from Theorem 5.1 to Theorem 5.3 is rather similar to 
that from Theorem 4.1 to Theorem 4.3. We note that the estimates done in 
the course of the proof of Theorem 5.1, with n > 3 there, can be summed up 
in the inequality 

|| x -(-+i)/^(^, x (-i)/2/)||^ < C(||/|| L oc)||/||^ , (5.40) 

where 

/ : = (f,X(l)A,X(f> + ,X(f)^) . 

The minor modifications of the proof of Theorem 5.1 needed to obtain (5.37) 
and the estimate (5.38) on (x + 2r)d T f are identical to the ones described in the 
proof of Theorem 4.3. The estimate on ||a;5 T + ||x,oo(n t j is obtained directly 
from Equation (4.49) and from (5.40). The estimate on ||a;d T (/>A||i,oo(n j is 
obtained from the (4.38a)-equivalent of the first of Equations (4.14). Next, for 
k > n/2 + 2 Equations (4.42) and (5.40) give 

e + {4>-) G ^f-i C , (5.41) 
Differentiating Equation (5.26) with respect to x gives 

d x 4>_{T,x)=d x (/>-(0,x + 2T)+ f \d x e + {<t>-)){a,2{T-a)+x)da , (5.42) 

J o 

which together with (5.41) implies, by straightforward integration, 

x(x + 2 T )\d x (f)-(T,x)\ < C (5.43) 
This, (5.41), and the identity 

d T (t>- = (d T - 28 x + 2d x )<f>. = e+(0_) + 2d x <f>- 
establish (5.39). □ 
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5.2 Estimates on the time derivatives of the solutions, n > 3 

To control the time derivatives of the solutions, as in Section 4.2 we introduce 
an index m which counts the number of corner conditions which are eventually 
satisfied by the initial data at the "corner" r = x = 0. As before we make a 
formal statement only for solutions of the wave-map equation (5.1), it should be 
clear from the proof that an analogous statement holds for solutions of (4.38)- 
(4.39) under appropriate conditions on the coefficients there. 

Theorem 5.4 In dimension n > 3 let N B m > 0. Consider a solution / : 
^x ,r t —>■ R of Equation (5.1) satisfying 

r 

||xe + Cr)|| L ^ ( a^ T j + \\ xX if a h°°(n X0 , T ,) 

i=i 

+ ll/ a |U~(n^) + Hz^/Iloo^) < oo , (5.44) 

and suppose that 

0<i<m + l d\Jy xofi G ^ ro +i_ i (E X0)0 ) , (5.45) 

0<i<m «9a4/°k ,o G J^ +m -i(^xo,o) , (5-46) 

with some > § + 2, -1 < a < -1/2. Then for < r < r* and for < i < m, 
we have 

0<j + i<A; + m — n/2 

[(r + 2x)5 r F47 a | Sxo , r G ^(S^,,) n ^& m+ i-i-i( S *o,r) , (5.47a) 
0<j + i<fc + m- n/2 - 1 

a x [(r + 2x)d T ydlf a \^ T G ^fc+ m _i_j(S X0)T ) , (5.47b) 

and 

< p < k - n/2 [(r + 2x)a T ]^r +1 7 a |E ;co , T € ^_ P (S X0)T ) , (5.48) 
with r-independent bounds on the norms. 

Proof: The proof is an inductive application of Theorem 5.3, as in the proof 
of Theorem 4.4, and will be omitted. □ 

5.3 Estimates on the time derivatives, n = 2 

In space-dimension two the following equivalent of Theorem 5.3 holds: 

Theorem 5.5 Let n = 2, consider the system (4.38)-(4.39), suppose that (5.31) 
holds for some constant C, with the nonlinearity G in Equation (4.38a) of the form 

G = x- z l 2 H{x^x 1 ' 2 ], x l ' 2 (j) A , x l/2 <p.,x z/2 <P + ) , (5.49) 

with G e _(^ A ) = (cf. Equation (4.41)), and with H having a uniform zero of order 
£ > 3 in the sense of (A. 30). Suppose that the initial data satisfy 
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/ a k ,o=^/ a k ,o G K a +1 n^ )(S XO) o), (5.50) 

Sx(n-2/°)|E 1S0 , e ^(S X0)0 ), (5.51) 

9r(n-^/°)|E <so ,o e (jr fc Q nL-)(s XOi0 ). (5.52) 

for some fc > 2, -1 < a < -1/2. Then: 



1. There exists t + > 0, depending only upon the constant C in (5.31) and 
a bound on the norms of the initial data in the spaces appearing in Equa- 
tions (5.50)-(5.52), and a solution f a of Equations (4.38)-(4.39), defined on 
a set containing £l X(hT+ , satisfying the given initial conditions, such that 

ll/ lk(nW < 00 • ( 5 - 53 ) 

2. Further, for any such that f a exists on fi X0)T „ with (5.36) holding with 
r+ = r*, we have for all < r < r* 

/ a k , T e « n M> k \i) (e xo ,t) , d r / a k , T g (jr£ n l°°) (£ X0 , T ) , 

with bounds uniform in r. This implies 

||xd r <£ + || L oo (QxojTt) + ||5 T 7 a || L oc (Qxo >rt) < oo . (5.54) 
If k > 4 then we also have 

|| (x + 2r)a r ^||L=o ( n <B0 , r j < oo • (5-55) 
If k > 4 and if d*f\v XOiT G then it further holds that 

||(x + 2T)S T ^_|| £ oo (n<e0iT ,) < oo . (5.56) 

Proof: The proof of point 1. is essentially the same as that of Theorem 5.1, 
with the modifications discussed in the proof of Theorem 4.3. We note that the 
key estimates (5.17) and (5.20) hold in exactly the same form here, similarly 
for Equations (5.29)-(5.30). The estimate on d T f in (5.54) follows from the 
definition of the norm in (5.53). The estimate on ||x<9 T <^ + ||£oo(ft xo t j is obtained 
directly from Equation (4.49) and from (5.17). To obtain (5.55) one needs 
to prove a bound on 8a4>-- This is obtained by differentiating (5.26) with 
respect to v A and using the already known uniform bound for G in Jtfj?, so 
that d A G e .^_ x C <*f Q . Finally, 

e+((x + 2t)3 t ^) = {x + 2t)8 t (e+(0_)) , 

and integrating as in (5.26) one finds 

(x + 2T)d T 4>^(T,x) = (x + 2r)<9 T <M0,x + 2t) 

+ / {(2(r - a) + x)8 T (e+(0_)) {a, 2(r - a) + x)} da . (5.57) 

JO 
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The term at the right-hand-side of the first line of (5.57) is bounded because of 
the hypothesis on d^f. Expressing e+(</>_) by the right-hand-side of (4.42), one 
immediately finds that all the linear terms that arise after differentiation with 
respect to r are in J^^-i or better, and therefore give a finite contribution when 
integrated upon. The contribution from the non-linearity G can be rewritten 
as 

X- 1 [Hjd T f + H^Or^ + H^drf- + H^xd^} , 

with appropriate functions which, by Lemma 4.6, all have a uniform zero 
of order I — 1 > 2 in their arguments. This easily implies that the coefficients 
(including the x~ x factor) in front of the r derivatives are in L°°, and since 
each of the r-derivative terms is in or better, the whole term is in 

c Co" -1 - This is sufficient to lead to a finite contribution in (5.57), 
and (5.56) follows. □ 

We finally arrive at the two-dimensional equivalent of Theorem 5.4; com- 
ments identical to those made in Section 5.2 apply here. The main difference is 
that in dimension 2 we need the L°° bound on d T f to obtain existence, which 
leads to the compatibility condition (5.62) on the second r derivatives of / when 
one attempts to iteratively apply Theorem 5.5. The proof is again identical to 
that of Theorem 4.4 and will be omitted. Let us just mention that one easily 
checks that the conditions spelled out below guarantee that the initial data for 
the inductive system of equations are in the right spaces for the iterative ap- 
plication of Theorem 5.5. Further, Equations (5.54)-(5.56) provide the a-priori 
bounds which guarantee that the existence time of the solution will not shrink 
at each iteration step. 

Theorem 5.6 In space-dimension two let N B m > 0. Consider a solution / : 
^x ,t, —> R °f Equation (5.1) satisfying 

H/l^(n W <°o. (5-58) 

and suppose that 

0<*<m di/°k , G (^Wi-in*f)(S XOl0 ), (5.59) 

0<*<m ^4/ | Sxoj0 G Jt%. m -i(2 xo ,o) , (5.60) 

dT +2 f a K,o G ^(Sxo.o), (5.62) 

with some k > 4, — 1 < a < —1/2. Then for < r < r* and for < i < m, we 
have 

0<j + i<k + m — 3 

[(r + 2x)d T Ydir\^ T G (^ Q +m+1 _ i _ j n 0(5^) , (5.63a) 
0<j + i<k + m — 3 

d x [(T + 2x)d T YdlJ a \^ T G ^U^-CExcr) , (5.63b) 
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and 

< p < k - 3 [(r + 2x)a r ]^- +1 7 a | S;co , T G ^_ P (S X0)T ) , (5.64) 
with r-independent bounds on the norms. 



□ 



5.4 Polyhomogeneous solutions 

We are finally ready to prove polyhomogeneity at J? of solutions of the wave 
map equation: 

Theorem 5.7 Let 5 = 1 in odd space dimensions, and let 5 = 1/2 in even space 
dimensions. Consider Equation (5.1) on W 1 ' 1 , n > 2, with initial data 

<9r/ a |{r=o} € n L°°)(M Xo ) , z = 0,1, n = 2, (5.65) 
7 a | {r=0 }G(^nL 00 )(M :Co ), d r / a | {r=0} G ^,(M I0 ) , n>3. (5.66) 

Then: 

1. There exists r + > such that / a exists on fi X0)T+ , with 

ll/ a |k°(n W < oo , n = 2, (5.67a) 

r 

\\xe + {J a )\\ L oo (nxQ T+) + ll^i/ a ||L«»(n B() , r+ ) 

i=i 

+ ||/ a ||L~(^ , r+) + H^(^ , T+ ) < oo , n > 3 . (5.67b) 

2. If the initial data are compatible polyhomogeneous in the sense that 

Vi G N 4/°(0) G L°°(M X0 ) , 

then the solution is polyhomogeneous on each neighbourhood fi X0)Ta , of 
on which / exists and satisfies (5.67) with r + replaced with r*. 

Proof: Existence of solutions follows from Theorem 5.1. Theorems 5.4 and 
5.6 give the time-derivative estimates which are necessary in Theorem 4.9. In 
order to apply that last theorem, we set 

*--(«)■ (5 ' 68) 

and 

v>i = cn, ^ 2 = w>-)- (5.69) 

Equation (5.2) takes then the form (4.58). As in Theorem 4.10, for n > 4 even 
we take <5 = 1/2, p = n + 3, q = n — 1; while for n > 3 odd we take 5 = 1, 
p = 2±3 > q = n^k. For n = 2 we set 5 = 1/2, p = 3, q = 1. The non-linearity 
here has a uniform zero of order 3, which is compatible with the hypotheses of 
Theorem 4.9, and the result follows by that last theorem. □ 

Acknowledgements: We are grateful to Helmut Friedrich for many useful 
comments on a previous version of this paper. 
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A Function spaces, embeddings, inequalities 



Throughout this paper the letter C denotes a constant the exact value of which 
is irrelevant for the problem at hand, and which may vary from line to line. 
Let M be a smooth manifold such that 

M = M U dM 

is a compact manifold 10 with smooth boundary dM. Throughout this work the 
symbol x stands for a smooth defining function for dM, i.e., a smooth function 
onM such that {x = 0} = dM, with dx nowhere vanishing on dM. It follows 
that there exists xq > and a compact neighborhood K of dM on which x 
can be used as a coordinate, with K being diffeomorphic to [0, xq] x dM. For 
< x\ < X2 < xo we set 

M X1 = {p <E M | < x(p) < xi} , (A.la) 
M XuX2 = {p € M | Xl < x(p) < x 2 } , (Alb) 
d~M Xl = {p£ M \ x(p) = xi} k, dM . (A.lc) 

In what follows the symbol f2 will generally denote one of the sets M, M XI , 
or M xljX2 . Any subset of M xo can be locally coordinatized by coordinates 
y % = (x, v A ), where the v A, s can be thought of as local coordinates on dM. We 
cover dM by a finite number of coordinate charts Oi so that the sets Q{ , where 

fii := (0,x ) x d , 

cover M XQ . We use the usual multi- index notation for partial derivatives: for 
= {Pi, . . . , 0n) £ N n we set d 13 = d{ 1 . . . dt ■ We will write d% for derivatives 
of the form d% 2 ■ ■ ■ dn" , which do not involve the x 1 = x variable. 

If & is an open set, for k G N U oo we let Cu{@) denote the usual space of 
fc-times differentiable functions on ff; the symbol Ck{&) is used to denote the 
set of those functions in C\ t [C') the derivatives of which, up to order k, extend 
by continuity to G. We emphasize that no uniformity is assumed in Ck(G), so 
that functions there could grow without bound when approaching the boundary 
of G. Nevertheless, the symbol || • \\c k will denote the usual supremum norm 
of / and its derivatives up to order k. The symbol Ck+\{G) denotes the space 
of fc-times continuously differentiable functions on G, with A-H61der continuous 
fe'th derivatives. 

For a e R, k G N and A € (0, 1], we define ^(Qi) (respectively ^ +x (^i) , 
^^(fij) , < ^ +A (^j)) as the spaces of appropriately differentiable functions such 
that the respective norms 

ll/ll«8»(no = sup \x~ a f(p)\ , 
pecii 

n-ii I, ,„ x(v)- a - x \f(.y)-f(.y')\ 

11/11^(00 = WfWv^) + sup sup — — - x , 

10 We use the convention in which manifolds are always open sets. Thus, a manifold with 
boundary does not contain its boundary as a point set. 
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= E ll^^/ll^nO) 

0<\(3\<k 

= Il/H^^) + E ll^^/ll^^) , (A.2) 

|/3|=fe 

are finite. Let & be an open subset of M, or a submanifold with boundary in 
M; for such sets we define: 



<e«{C) = SU P ll/lk«(Qin^) + \\j\\c k {GM xo/2 n0) > 

= su Pll/lk" +A (n<n^) + ll/llc fc+A (CM Io/2 n<?) • ( A - 3 ) 

We note that / G < ^ f A CT (^) if and only if G ^f +A (ft). 

We define the spaces J^ a (fij) as the spaces of those functions in i^ oc (f2j) 
for which the norms || • H^gy*^) are finite, where 

2 *?W= E / (A.4) 

0<\l3\<k JQi 

Here dz^ is a measure on <9M arising from some smooth Riemannian metric on 
dM. This is equivalent to 

£ / {x-^xd^dlff^dv, (A.5) 
0</3i + |/3|<fc 4 

and it will sometimes be convenient to use (A.5) as the definition of H/H^am.v 
For ^"s such that Qi C 6 the spaces ^^{G) are defined as the spaces of those 
functions in H l g c (ff) for which the norm squared 

Wff^(ff) = E WfWje? (fii) + ll/llff fc (tfnCM xo/2 ) ( A - 6 ) 

is finite. We note the equivalence of norms, 

l^o(^) ~ ll/ll ^- 1 ' 2 i t ff\ ' 



and that J^ x (M XltX2 ) = Hk(M Xl X2 ) for all a and A; whenever x\ > 0, the norms 
being equivalent, with the constants involved depending upon x\ and X2, and 
degenerating in general when x\ tends to zero. 

It is often awkward to work with coordinate charts, in order to avoid that 
one can proceed as follows: Choose a fixed smooth complete Riemannian metric 
b on M. Let x be any smooth defining function for dM, we let X\ be the 
gradient of x with respect to the metric b; rescaling b by a smooth function if 
necessary we may without loss of generality assume that X\ has length one in 
the metric b in a neighbourhood of dM. As before we cover dM by a finite 
number of coordinate charts Oi with associated coordinates v A ; the v A, s are 
then propagated to a neighbourhood of dM by requiring 

X 1 (v A ) = 0. 
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This leads to a covering of M XQ of the kind already used, and one easily checks 
that 

X 1 = d x 



in the resulting local coordinates. This gives then a globally defined vector d x 
on M XQ . For i = 2, . . . , r we let Xi be any smooth vector fields on dM satisfying 
the condition that at any p G dM the linear combinations of the Xi exhaust 
the tangent space T p dM. (If dM is a sphere S" 1 " 1 , a convenient choice is the 
collection of all Killing vectors of (S n ~ 1 ,h), where h is the unit round metric 
on S n ~ 1 .) Over the domain of a chart (v ) of dM, one thus has 



9a = £/A(v B )*i, (A7a) 

i=2 
n 

X t = J2 X f( vB ) d A, (A-7b) 

A=2 

for some locally defined smooth functions f\, Xf; clearly things can be arranged 
so that those functions are bounded, together with all their partial derivatives. 
We propagate the JQ's to M xo by requiring 

[X 1 ,X t } = 0, 

equivalently 

d x Xf = . (A.8) 

It follows that (A. 7) still holds with x-independent functions. For any multi- 
index (3 = (Pi, fe, ■ ■ ■ , Pr) G N r we set, on M Xo , 

®Pf = X{ X X^ ■ ■ ■ X&f = d^X^ ■ ■ ■ X&f . (A.9) 

It follows that we have 

i^(^ )~ E w^^fWvw**) ' 

0<\8\<k 



0<\f3\<k 



V / (x- a+ *@Pf) 2 — du 



(where denotes the fact that the norms are equivalent), etc. Here, \(3\ = 
ft + ... + &. 

There is a useful way of rewriting || • \\j?*/ Mx \ which proceeds as follows: 
for / G J^ k a (M X0 ), s G (1, 2), and n G N we set 

fn(s, v) = fix = x ^,v) ; (A.10) 

letting rs denote again equivalence of norms one then has, after a change of 
variables, 

Wo) = E E / , , „ \*- a+ ^f(*,v)\^dv 
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£ E ^ / \^fn(s,v)\ 2 dsdu 
n>10<|/3|<fc Ah2]xdM 

^ 2 2na \\fn\\H k ([l,2]xdM) ■ ( A -H) 



— 2a r,2na || f \\2 

n>l 



More precisely, we write A ~ i3 if there exist constants Ci , C2 > such that 
C\A < B < C2A. In (All) the relevant constants depend only upon a and k. 

It turns out to be useful to have a formula similar to (A. 11) for functions in 
J^(M X2jXl ); this can be done for any x\ and X2, but in order to obtain uniform 
control of certain constants it is convenient to require 2x2 < x\. For such values 
of x\ and X2 we let n§{x\,X2) G N be such that 2n ^+i < ^2 < For n € N, 
n > 1, and for any / : M X2jXl -> we then define /„ : (1, 2) x <9M -> by 

n<n , f n (s,v) = f{xi — ,v) , 
n = n + l, f n {s,v) = f(x 2 s,v) , 

n>n + l, f n = 0. (A.12) 

(This coincides with the definition already given for M Xl , when this set is 
thought of as being an "M X2)Xl with X2 = 0", if we set no = +00.) A cal- 
culation as in (A. 11) shows that for any 2x2 < x i < x o, there exist constants 
C\ and ci, independent of xq, x\ and X2, such that for all / G J^ a (M X2;X1 ), 

ClX r 2 " E{ 2na H/™llH fc ([l,2]x9M)} 2 

n 

^ H/HW^) ^ ^l^ 2Q E„{2 na ||/n|k fc ([l, 2 ]x 9 M)} 2 • (A.13) 



Equation (A. 11) leads one to introduce 11 spaces that arise naturally from 

weighted Sobolev embeddings, cf. Equation (A. 25) below: we define 

\\f\\%l +x {M X0 ) = x o 2a ^2 22na \\fn\\c k+x {[i,2]xdM) > ( A - 14 ) 

n>l 

/ n as in (A. 10), and we set 

^k+x(M X0 ) = {/ € WM X0 ) I ||/||^ +a( m xo ) < 00} . 

Clearly 

^k + x(M X0 ) C <^ +A (M Xo ) . 

Since the general term /at, as well as sums of the form T> n >Nfn, of a convergent 
series tend to zero as N tends to infinity, for / £ S^ +A (M :E() ) we actually have 

^11/11^(^=0. (A.15) 
We have the trivial inclusion, 

0^i) c ^( M ^)- (A. 16) 

lx The symbol ^ might suggest to the reader that we specifically have Besov spaces in mind; 
this is not the case, and we hope that the notation will not lead to confusion. 
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The fact that the inequality a' > a in (A. 16) is strict has various annoying 
consequences, which are best avoided by introducing yet another space — the 
space ^ of functions in H^ oc (M X{) ) for which the norm squared 



V / \x- a+ ^$>Pf(x,v)\ 2 — du 

(A.17) 

is finite. We note that ||/||^«(m xo ) is equivalent to 

x^sup {2 na \\f n \\ Hk{[1 , 2]x9M) } , (A.18) 

n>l 

with f n (s,v) = f(^£,v), as in (A. 10). To define the ^*(M X2)X1 )'s, assuming 
again that x 2 < xi/2, we let 7 n (xi,x 2 ) be defined as 

n<n , I n = (2- n Xi,2 1 - n xi) , 
n = n + 1 , 7 no+ i = (x 2 , 2x 2 ) , 

n > n + 1 , = , (A.19) 

where no is as in (A.12). For all / G H k oc (M X2 , Xl ) we set 

W^^E £ /„ , „ y a+ ^f? d idv) (A.20) 



(we identify (a, 6) x <9M and M a ^). Similarly to (A. 13), there exist constants 
c 2 and C 2 , which do no£ depend upon xo, x±, and x 2 , such that for all 2x 2 < 

Xl < x , 

C 2 xf°SUp ||/n||ff fc ([l,2]x9M) < ll/II^M*,*,) < SU P II /« II H fc ([l,2] xdM) ■ 

n n 

(A.21) 

We have the obvious inequality 



|s?«(n) < ||/||^«(Q) , (A.22) 

together with the modified version of (A. 16), 

a >> a ^ A c^i (A.23) 

in particular the function (x,f) — > x a is in & k (M XQ ). 

If 5fc denotes a space of functions, where k G N is a differentiability index, 
we set 

Soo = n^^s^ , 

e. 5 ., «S = n km ^, etc. 
We note the following: 

Proposition A.l Let fl = M, or f! = M Xl , < xi < xo, or = M X2)Xl , 
2x 2 < xi < x , and let Jf k a = Jf k a (n), etc. For k' G N, A G [0, 1], < fc' + A < 
/c — n/2 G"NorO</c' + A<A; — n/2 G N we have the continuous embeddings 

C 38% +x C , C ^ C , (A.24) 
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and there exists an ^-independent constant C such that we have 

V/ G ^ ll/ll^ +A (n) < C||/||^- (n) , (A.25) 
V/ G ^ H/H* fc V A (n) < CH/b fc «(n) • (A-26) 

Proof: (A.25)-(A.26) follow immediately from (A. 11) and (A. 13), together 
with the standard Sobolev embedding; the remaining inclusions in (A. 24) are 
trivial. □ 



All other inequalities involving Sobolev spaces have their counterpart in 
the weighted setting; we shall in particular need various weighted versions of 
the Moser inequalities. The reader should note the different weights for the 
members of Equation (A. 31) below — this shift of weights in this inequality is 
the key to our handling of nonlinear equations. 

Proposition A. 2 Let $7 = M, or Q = M Xl , < x\ < xq, or Q. = M X2jXl , 

2x 2 <x x < xq, and let J? k a = J^ k a (ty, etc. 

1. There exists a constant C = C(a, a', (3, k, xi) such that, for all / G J^'ntftf 
and geJtfn ^ +f3 ~ a ' , we have 

\\fg\\^<c{\ 

Further, V | 1 < k, 



\9\\je£ 



+ 



^«'llfll^a+/3-a') 



(A.27) 



\x-n&>(fg) - (x^^f)g\\^ < C \\\f\\v °\\g\\^ + 

( \\ X 9xg\\^a+/3- a ' + ||^g||^a+/3-a' j j J (A.28) 



where the vector fields X are defined in Equation (A. 7). 



2. Let F G C k {M x R N ) be a function such that for all B G R + there exists a 
constant C\ = C\{B) so that, for all p G R N , \p\ < B, we have 



\F(; 



< Ci . 



Then for all a < 0, j3 G R, and B G R + there exists a constant CiiB-> k, a, j3, x\) 
such that for all R^-valued functions / G J^ a_/3 (J2) with Ha^/Hx^n) < B 
we have 



FCx' 3 /) <C 2 (1 + H/ll 



(A.29) 



Further, if F has a uniform zero of order Z > at p = 0, in the sense that for all 
B G R there exists a constant (7(1?) such that for all \p\ < B and < i < min(/c, I), 



#F(;p) 



dp 1 



«J-i(^-o) 



(A.30) 
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then for all a € K, (3 > 0, there exists a constant Cs(C,l, k,a, (3, B) such that, 
for all / G ^"^(fi) with ||/|| L oc (n) < B, we have 



(A.31) 



Remark: The hypothesis (A. 30) will hold if F is e.g. a polynomial in p with 
coefficients of p> vanishing for j < I, and being functions belonging to ff® for 
3>l- 

Proof: We shall give a detailed proof of (A. 29) and (A.31), the inequalities 
(A.27)-(A.28) follow by an analogous argument using [35, Volume III, p. 10, 
Equations (3.21)-(3.22)], cf. the calculation of Proposition A. 3 below. Let, 
similarly to (A. 10), 

t? I \ T?ll X S \ ( X Q S \0t, X S \\ 

F n (s,v) = F y(x = —,v); (— ff(x = —,v)) ; 
from Equation (A. 11) we have 

\\F{-, xl3 f)\\ 2 Jif k a (M X0 ) ~ x o a X/ 22na H^"ll^fc([i,2]xaM) • 



(A.32) 



n>l 



We have the obvious bound 



sup 

[1,2] xdM 



Further the partial derivatives of (s,v) — > F n (s,v,p) with respect to s and t> at 
p G 1^ fixed, |p| < M, can be bounded by a constant depending only upon 

sup ||F(-,p)|U 0(M } . 

|p|<M fc 

The usual Moser inequalities [35] [Volume III, p. 11, Equation (3.30)] give 

H-^ll!r fc ([l,2]xdAf) - C (i + 2 2n/3 H/™ll/f fe ([l,2]x9M)) ' 

with / n as in (A. 10), and with a constant C depending upon A; and M. Inserting 
this in (A.32) one obtains (recall that a < 0) 



\ F ('^ xl3 f)\\ 2 ^(M X() ) - C^2 2na (l + 2 2nl3 \\fn\\ 2 H k ([i t 2]xdM)) 



n>l 

< c(i + 



(A.33) 



This establishes (A. 29) for = M xo , and (A.29) with Q, = M readily follows. 
The remaining fi's are handled in a similar way. 
To establish (A.31), we note the inequality 



d^ +i F n (;p) 
dy^dp 1 



<C\p\ 



max(Z— i,0) 
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which follows from (A. 30) when \ j\ + i < k. Letting y stand for (s, v) G [1, 2] x 
dM, it then follows that for \a\ < k we have 



|S"F„| = 



|7l+kl|+-+|<Ti|=|cr| 



X ^# aCT1(S / " ) ---^ (S /n) 

< 2-'^C \d ai (^fn)\---\d' 7i (sPf n )\. 

ki|+-+ki|<|<r| 

The usual inequalities [35, Volume III, Chapter 13, Section 3] give 

\\Fn\\H k ([l,2]xdM) < C(fe, ^) 2 ~' /3n ||/n||// fe ([l,2]x9M) i 

for some constant C(k, M), and one concludes from (A. 32), as in (A. 33). □ 

We have the following sharper version of (A.27)-(A.28): 

Proposition A. 3 Let Q = M, or $7 = M Xl , < x\ < xq, or = M X2jXl , 
2^2 < x\ < xq, and let = Jif^(Q), etc. There exists a constant C s = 
C s (a, (3, k) such that, for all / e ^ a n ^ and <? G Sff n we have 

\\f 9 \\^ < C s (\\f\y s \\g\\^ + H/ll^a \\g\\ % e) , (A.34) 
Moreover it also holds that 



V| 7 |<fe, \\x^^(fg)-(x^^f)g\\^ a 



< C \\f\y«\\g\Le + H/H^ ||x^|| , + ||A^|| p (A.35) 





i=2 



where the vector fields X are defined in Equation (A. 7). 

Remark: A useful, though less elegant, inequality related to (A.34) is 

V| 7 +<r| < k \\x^(@y)x^(@°g)\\^ + e < C s (\\f\\^ \\g\\^+\\f\U k " Nl^) ■ 

(A.36) 

Proof: We will prove (A.35), the proof of (A.34) is essentially identical. When 

Q = M xo we do the rescaling f n (s,v) = 9n(s,v) = g(®&,v), we then 

have, for all \ j\ < k, 

\\x^^(fg)-{x^^f)g\\ 2 ^ +[) 

« H (a+P) E 22n(a+/3) ll^(/^) " ^fn)gn\\ 2 Ho{[1 , 2]xdM) 



< Cx~ 2 ^ 



E 2Ma+/3) (ll/nlli-ll^ll^ + \\fnt Hk J\®9n\\l~) 
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(A.37) 



(In the third line above we have used the inequality [35, Volume III, p. 10, 
Equation (3.22)].) The case 0, = M follows immediately from the above; the 
case £1 = M X2X1 is treated similarly using (A.12)-(A.13) and (A.19)-(A.21). □ 

Similar results can be proved in weighted Holder spaces: 

Lemma A. 4 Let $7 = M, or ft = M Xl , < x\ < xq, or £1 = M X2jXl , 2x2 < 
xi < s . and let fc£ = Let / G n <*f<f and fl G %? fc 7 n ^ with 

a + 5 = 7 + /? = <T. Then we have fg G ^ and 

ll/ffll^ < Qdl/II^I^H^ + Nl^||/lk«) , (A.38) 

Proof: The proof is very similar to that of Propositions A. 2 and A. 3. We use 
the same conventions as in (A. 12), (A. 19). We have ~ sup ra 2 nr7 \\f n g n \\ Ck ^, 

where 

u = [1, 2] x dM , (A.39) 
similarly for / and g. The interpolation inequality [27, Appendix A] 

Wfn9n\\c k {uj) < C(ll/n||oo||ffn||c fc (w) + llffn ||oo ||/n||c fc (u;)) 

leads to the conclusion. □ 

We have the following ^ equivalent of the second part of Proposition A. 2, 
with a similar proof, based on Lemma A. 4: 

Lemma A. 5 Let F be a function satisfying the hypotheses of point 2 of Proposi- 
tion A. 2, with a uniform zero of order I'm p in the sense of Equation (A. 30). Then, 

for any e > 0, p G R and / G n L°° we have F(.,x e f) G and there 

exists a constant C depending upon ||/||l°° such that 

\\F(.,x*f)\\^ + u<C(\\f\U\\f\\^. (A.40) 

The space of polyhomogeneous functions £^ p hg = ^phg(M) is defined as the 
set of smooth functions on M which have an asymptotic expansion of the form 

oo Ni 

/~EE/^"' I"'-'-. (A.41) 

i=0 j=0 
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for some sequences rii,Ni, with rij /* oo. The poly homogeneous expansions of 
the introduction are of this form if r there is replaced by 1/x; this corresponds 
to the conformal transformation of Section 2, which brings "null infinity" to 
a finite distance. We emphasize that we allow non-integer values of the jij's; 
however, we shall mostly be interested in rational ones, as those arise naturally 
in the problem at hand. Here the symbol ~ stands for "being asymptotic to": 
if the right-hand-side is truncated at some finite i, the remainder term falls off 
appropriately faster. Further, the functions are supposed to be smooth on 
M, and the asymptotic expansions should be preserved under differentiation. It 
is easily checked that the space &? p hg is independent of the choice of the function 
x, within the class of defining functions for dM. 

B ODE's in weighted spaces 

In our handling of PDE's below we will need ODE estimates to obtain informa- 
tion about solutions, we thus begin with some a priori estimates in weighted 
spaces for ODE's. While the results are well-known in principle, and easy to 
prove, we present them in detail here because their precise form is necessary 
for our arguments later in this work. For a vector w we denote by or by 
\w\ the usual Euclidean norm, while for a matrix b the symbol ||6|| denotes its 
matrix norm. 

B.l Solutions of d T (p + b(p = c in weighted spaces 

Let O be an open subset of dM, which might be the whole of dM, or a coordi- 
nate patch of dM with coordinates v A , whichever appropriate in the context; 
we set 

U X2 , Xl = (x 2 , Xl )xOx[0,T], (B.l) 
y x2 , xi = (x 2 ,xi) x O , (B.2) 

with < X2 < x±. The time variable r will usually be the last variable, so r 
will run from [0, T] whenever U X2;X1 is involved. Strictly speaking, U X2;X1 should 
carry an extra T index, but we have not done that in order not to overburden 
notation. To avoid ambiguities we emphasize that the spaces c ^ , ^{Ux 2 ,x\) in the 
Proposition below are defined as in the previous section, with the v A variables 
there corresponding here to some local coordinates on O together with the time 
variable r; the time derivative d T should be understood as a one-sided one at 
r = and at r = T. 

Proposition B.l Let a € R, b G ^(U X2tXl ,End(R N )), c G tfg(U X2tXl ,R N ), 
then the unique solution tp of the equation 

8 T ip + bip = c , (B.3) 

with initial data <p = <p\ T=0 £ tf£(y x2iX1 ,R N ) is in tf£(U X2tXl ,WL N ) with 

M^(U X2 , X1 ) < C (n,N,k,T, Xl ,\\byo {Ux2 ^ (\\(f>\W^ 2 , xl ) + \\4^(U X2 , X1 )) ■ 

(B.4) 
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We also have the estimates 

MT)\\ Wx2tXl) < Ce^ (M0)|| WW +^ T e- ||6||ooS ||c( S )|k« ( ^ 2 ^) ds 

(B.5) 

II^IIWW ^ C e c " b »^x (11^(0)11^^^) +^%- c l 6 l^||c( s )||^ ( ^ 2 ^ i) d S 

+ / T e (1 - c)|fc| - l&OOIk^w (mo)!!^^) 



+ ^ S e-l b l- t ||c(t)||^ ( ^ i) ^) da) . (B.6) 



Remarks: 

1. Analogous results in spaces can be proved by similar arguments. 

2. An a-priori estimate in weighted Sobolev spaces for (B.3) follows from 
Proposition 3.1 below by setting E^d^ = d T ® id and L = ip = b = there. 

Proof: Let k e W, and let j3 = /3 2 , • • • /?n) be a multi-index with < k; 
d 13 ^ verifies the equation 

d^ip = -d p {b^)+d^c. (B.7) 

Let e > and set 

\ V2 

e(.,t,e) = [e+Yl ^ 2(/3l ~ a) (dV &*<p) 
\ \P\<k 

E(t,e) = \\e(.,t,e)\\ L oo { ^ X2Xi) . 
When k = one easily finds 

<9 r e < + |c| , 
and (B.5) readily follows. For A; > we have 

d T e = -Y, x m - a) (d T d^, &<p) , 
e \/3\<k 

- x 2 ^- a ^(-bip + c)\ \d^if \ , 
6 \/3\<k 
C(k,n) 



< 

6 \a i. 



< 



e 

< C(k,n)(\\bip\\^ { ^ 2 xi) + \\c\W s[ ^ X2 Xi) ) , 

where C(k,n) is a constant depending upon k and the space dimension n, and 
which arises from the inequality Yli=i \ a i\ — \ft> vSl^F f° r an y rea ^ sequence 
(dj). The weighted interpolation inequalities, Lemma A. 4, imply 
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where C is a constant which depends upon k, N and n. It follows that 

d T e < C (\\b\\ L o°^ X2 Xi) \\v\W S (^ xl) + \H^ k {.y X2tX1 )\W\Wz(.?' X2 , Xl ) + l|c|k« ( ^ 2 xi) ) 

< C (ll&Hoo^M) + l|b||^o ( ^ 2iXi) ||^||^a ( ^ 2iXi) + \\c\W s( ^ xi) ) , 
with perhaps a different constant C. By integration we obtain 

e(t) < e(0) + C f (||&|| 0O £7( s ,e) + IIK*)!!^^^^) Il^s)!^-^^^) + llcWH^^^^)) ds , 
from which we deduce 

£(t,e) < £7(0, e) + C f (||6|UE( S ,e) + ll&^ll^^^)!!^)!!^^^) + IK^Ik-pw) <fe 



Using Gronwall's Lemma and letting e — > one obtains 
£7(t,0) < e c||fe|l °°*£;(0,0) 

The estimate (B.5) for WvW'gg'^^ ) inserted in the last inequality leads to 
Equation (B.6). The time-derivative estimates follow immediately from the 
above and from the equation satisfied by ip. □ 

B.2 Solutions of d x (f) + b(f) = c in weighted spaces 

All the results in this section, as well as in Section B.4 below, remain valid if 
we replace the set U X2;X1 defined in Equation (B.l) with ^ 2)X1 defined in (B.2) 
- the time dimension does not play a preferred role in the current problem. 
We start with the following elementary result; the point is to ensure that the 
relevant constants are x 2 independent: 

Lemma B.2 Let g G tf£(U X2 Xl ,R N ), < x 2 < x lt then / defined for a > -1 
by 

r-x 

f(x,v A ,r) = / g(s,v A ,r) ds 
is in V° +1 (U X2 , X1 ,R N ), with 

Wf\\^ +1 (u X2 , xl ) < max {^^Jj} \\9\\v?W*2.«i) ■ 
Similarly f 2 defined by 

f 2 (x,v,T) = - g(s,v,r) ds 

J X 

satisfies 

(1 + (Inxfr^h e ^(U X2 , X1 ) for a = -1 , 
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/2€ 



? min{a+l,0} 



{U X2jXl ) for a < and 



with 



II/2II ^mi^a+i.o}^ ^ ^ < max < 1, 



Proof: We have the trivial relations 



s a ds < 



1 



1 + a 



1 + a 



\9\\V{?(Ux2,x 1 ) 



1X2 



a + 1 



for a > — 1 , 



ds = In xi — In x , 



as well as the commutation rules: 

d x g dx 

J a 
i-x 

d v A / gdx 

J a 



d T I g dx 



g{x) , 

/•X 

/ d vA g dx , 

J a 

d T g dx . 



Note that 



l^/U^W™) + £ ll^/ll^^..,) , (B.8) 

0<i+|<5|<fc 



with \\d x fy«_ i{Ux2Xi) = \\g\\^_ l( u X2 , xl )- To estimate fyd^f one writes 
Mf\ < [ X \dld 5 v g\ ds , 

J XI 



1X2 

< [ x \\did 5 vg y (? s a ds , 

J X2 



< 



1 



a + 1 



The results for f 2 are established in a similar way. 
We shall use the following notation 

j X2 ={x = x 2 } , 



□ 



(B.9) 



with the range of the other variables being in principle clear from the con- 
text; this is the equivalent of the set dM X2 of Equation (A.l) when the set-up 
described there is assumed. 

Proposition B.3 Let < x 2 < x\, suppose that b G ^ e (U X2jX1 ,End(R N )), 
< e < 1, c € tfg(U X2jX1 ,R N ), and let be a solution in C l £ c '(U X2jX1 ) of the 
equation 

d x 4> + b^> = c. (B.10) 

Then the following hold: 
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1. If a < -1 , then <p G ^ + (U X2jX1 ) and we have, for a + 2 - e / and for 
^2 < ^3 < #i small enough so that C(||6||^.- e , £3) < 1, x 3 7^ 0, 

ll^k+Ws) - l-C(||6||^,X3) ( ^ 1|l ^ llco( ^) + ^^ll c ll^(^-3)) • 

(B.11) 

where 

Moreover, if £2 < X3 < xi is small enough so that C% C(||6||^- e , x 3 ) < 1, 
where C, is the constant in the interpolation inequality (A. 38), then 

H^k a+1 (Z4 2 ,* 3 ) ^ Ca(||6||^-e,Ci,X 3 ) (||0(^3)||c fe (^ 3 ) + l|c|k« (Wx2 ^ 3) 

+ ll 6 H* fc - e W.2,-3) (ll ^ a:3)llc °(^3) + ll C ll^(^ 2 ,^))) ' 

(B.13) 

with C a (\\b\\<g,-e , Ci, X3) an increasing function in the first and third variable. 

2. If a = 1, then (1 + (Inx) 2 )- 1 / 2 ^ G ^(U X2yXl ). 

3. If a > —1, then X2 = \im x ^ X2 4> is in Cfe(j^ 2 ), with 

4> - <j> X2 G <- e (Z4 2lX1 ) + K +1 (U*2,*i) > (B.14) 
0G^ Q+1 (Z4 2):ri ) if <^ 2 =0, and 

1 / \ 

l^^-W^- a) " l-C'(||6|| y -.,x 3 ) (j^^-a) + C H WW J 

(B.15) 

for X2 < X3 < xi small enough so that 

^(11*11^.^3) :=^\\b\\^ e{Ux2X3) <l. 
Moreover for x 3 small enough so that CiC'(\\b\\^,-e,xs) < 1 we also have 

11^11*2(^2,-3) ^ Cadl&ll^-^Ci^s) (ll^3)||c fe (^ 3 ) + l|c|k«(w* 2 ,* 3 ) 

+ ll 6 ll^r(^2^)(ll^ 3 )ll c o(>x 3 ) + ll c lk«(w* 2 ,* 3 ))) > 

(B.16) 

with an increasing function in its first and third argument. 

Remarks: 1. The inequalities above are standard when X2 > and when the 
constants are allowed to depend upon X2, regardless of whether or not X3 can 
be made small. As already mentioned, the point here is to make sure that the 
constants do not blow up as X2 gets small. 
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2. In case 2. log-weighted estimates are easily derived; they will, however, 
not be needed in what follows. 

Proof: 1. For simplicity, we will write ^ for ^(U X3tX2 ). Let <f> be a (local) 
solution of (B.10), corresponding to initial data at {x = x\} in Cfe(J^ 1 ). For 
a > set 

e a (x,v A ,r) :=(a+J2 ^ (^l^)) 1 / 2 , 

|/8|<* 

and e := eo- Let X3 G]x2, xi[n]0, 1] be such that \2+a-e\ W^\\v € < ^ ^ e ^ ave 
for all X2 < x < X3, 

-d x e a = --Y^I3ix 2 ^~ l {d^\d^) I 

-l^x 2 ^(^^) II, (B.17) 

6fl l/3|<fc 

Since /3i is non-negative we have — d x e a (x, v A , r) < II ; further 
II = — x 2ft (^(6^-c)|^) 

< — {\x Pl d p c\ + (^^(^l) \x^dP<t>\ 
6a \/3\<k 

< Yl \x Pl d p c\ + |a;^^(60)| . (B.18) 



Clearly 

J2\x^d^c\ = x a ^\x- a+ ^d f3 c\ 

< x a \\c\\^ , 

\x Pl d^(b<j))\ = x a+l ~ e Y \x- a - 1+t+Pl d p (b<j))\ 

\P\<k \&\<k 

which gives 

-d x e a < x a \\c\\^ . (B.19) 

k k 

Consider, first, the case k = 0; in this case (B.19) reads 

-d x e a < x a \\cy° +x a+1 -*\\b\\^ r \\cP\\ K+ i , 

which, after integrating over [x^, x] and passing to the limit a — > 0, gives (recall 
that a < — 1) 

e(x,v\r) < e(xs V,r)+(-^l y + ^)|| C ||^ 
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c a+2 ~ e x a+2 ~ e 



_ (2 + a-e) (2 + a-e) 
< Mco(^ 3 ) + ji^ll c K a 

(a+2-e ™a+2-e \ 
7——^ r - -r r . (B.20) 

(2 + a-e) (2 + a-e) J " "^0 imi *o v ; 

Suppose for the moment that a + 2 — e < 0; Equation (B.20) yields 

e(x,v A ,r) < U\\ Co ^ 3 ) + |T^ll c lk" + \ 2 + a-e\ M ^ ^'k +1 ' (R21) 
and since x _1 ~ a < x$ 1 ~ a < 1 we obtain 

x- a - 1 e(x,v A ,r) < x^UWco^) + ll c ll«? + \ 2 + a- e| W h K^K +1 

On the other hand, ifa + 2 — e>0 then 

x a+l x a+2 ~ € 

e(x,v A ,r) < ll0llco(^ 3 ) + |n^ll C ll^ + (2 + Q _ e) )ll b ll^^ll^ +1 ' 



which gives 



x a ^(z.vV) < V a |l^llco(^ 3 ) + ^^l|c||^ 



,1-e 
13 

|2 + a — e| " 



+ ^ rllft|l*-«ll0IU + i • 



The inequality H^H^a+i/^ ^ < sup^ ^j x a e shows that in all cases we 
have 

with the constant as in Equation (B.12). Consider, now, any < k E N; 
Equation (B.19) and the interpolation inequality (A. 38) give 

-d x e a < x a \\ c y« +x a+1 ^a(\\b\\^ r H\\^ + \\b\\ % -44>\\^) ■ 

An argument identical to the one before, considering separately the cases a + 
2 - e > or < 0, leads to 

•r' - 1 - GC(W„ r , I3 ) (^'""IMIw.,) + jrb^lMk 

+ i-GC(Sii y -^ ) |2+r.| l ' t| k-Wk« 



'0 



C(x 3 ) 



< ^ ^,„,m r ( Il0llr7,c^ 3 ) + \\c\\<tfg 



l-CiCqbWv-;**) 





C ? ; „, „ / _ Q _i r ,,, i 



1-C||&|U- E ,*3)" V d ' |l + a 



^3° Wco(^s) + 771- jll c lk Q 
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which gives (B.13). We have thus shown that (f> £ (U X2:X3 ); the property 
that (f> 6 c t?h +1 (Mx2,xi) immediately follows. 

2. The proof is identical, except for a few obvious modifications in the 
calculations. 

3. To obtain the L°° estimate, we start from (B.17)-(B.18) with k = 0, 
which upon integration and passing to the limit a — > gives 

x a+1 x 1 ~ e 
e(x,v A ,r) < e (x 3 ,v A ,T) + r^\\c\\^ + ^J\b\\^-4<pyo , 

from which we deduce 

\\<f>h°°(u X2 , X3 ) < II0IIl-(^ 3 ) + ^rlNk? + T37ll 6 H%- E 11^11^(^2,^) > 

and (B.15) follows. The proof of (B.16) is similar to that of the analogous 
statement in point 1. From what has been said it can be seen that (f> X2 = 
lim x ^ X2 (j) exists and is in Ck{J^ X2 )- It remains to show that <p — <p X2 satisfies 
(B.14). When b is a multiple of the identity, we can integrate (B.10) to obtain 

<f>(x,-) = X2 (-)e-^ b{s >- )ds + f ef* b ^ ds c(y,-)dy , (B.22) 

J X2 

from which the result easily follows. The general case can be established by 
manipulations similar to the previous ones. □ 



B.3 Poly homogeneous solutions of d T (p + b(p = c 

We pass now to an analysis of ODE's with poly homogeneous sources. The 
results here have an auxiliary character, and several of them are rather ele- 
mentary; they will be needed to handle the real problem at hand, with partial 
differential operators. Let O be an open subset of dM, we set 

U X1 =]0,xi] x O x [0,T] . (B.23) 

It will be seen in Sections 4 and 5 that 12 integer space-dimensions force us to 
consider poly homogeneous expansions with half-integer power of x; in order to 
account for that, we introduce an index 

-i 

where d is a non-zero integer, d £ N*. We will mostly be interested in the 
case d = 1/2 or d = 1, however other values are also possible in the formalism 
here. Results analogous to the ones below hold for the general poly homogeneous 
expansions of Equation (A. 41), which can be established by similar methods. 
We find it of interest that a consistent framework can be obtained in the setting 
considered below: 



12 This is due to occurrence of the factor fi' n in equations such as (4.2). 
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Proposition B.4 Let j3 6 R and consider the system 

d T ip + bip = c , (5.24a) 

V Ni 

(p\{ T=0 }(x,v) = <p(x,v) = ^2^2 x tS In 3 x <pij{x,v) + (p p5+ p +t {x, v) , 

(5.24b) 

€ Coo({r = 0}) , w,) +£ etf Wf ({T = 0}), (5.24c) 

with 

p ^ 

b(x,v,r) = x* 5 In 7 x &jj(x, t>, r) + & p< 5 +e (x, t>, r) , (5.25a) 

j=0 jr=0 

b p5+e e^ +£ (U xl ) , bij € , (5.25b) 

P AT" 

c(x,-y,r) = x^^^x^ln 5 x Cij(x,v,r) + c p< 5 +/3+(: (x, u, r) , (5.25c) 

i=0 jr=0 

E ^ +P+e (U Xl ) , Cij e Coo(U xl ) , (5.25d) 
where < e < (5, and (iVj), (AT'), (iV") are sequences with integer values, and with 

bGL°°(U Xl ) . 
Then the solution ip takes the form 

p M, 

99(x, u, r) = x^ ^ ^ x lS W x <pij(x, v, r) + ip pS+/3+e (x, v, r) , (B.26) 
i=0 j=0 

with ipij G Cqo^J, Mfc is an integer sequence and tp p $ + p +e e ^oo 

To prove the proposition we shall need the following lemma: 

Lemma B.5 Under the hypotheses of Proposition B.4, suppose that in addition 
we have 

<Pp8+/3+e = bp8+e = CpS+/3+E = . 

Then for any e g]0, 6[ we have 

V Mi 

^ = xP Yl Yl ^ lnj x + fpS+p+e > (B.27) 
i=0 j=0 

with ipij G Coo(U Xl ), (p p s+i3+e £ ^So +l3+e {U X1 ) , for some integer- valued sequence 
M fc . 
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Proof: Inserting (B.27) in the equation (5.24a) and tracking the coefficients 
in front of x tS In 7 x one finds the following set of equations: 

M = max{iVo, N%} , M i+1 = max{ max M k + N'^ k , N" N i+1 } , 

0<k<i 
i mm{N' k ,j} 

i G [0,p] , j e [0, Mi\ , drVij + ^2 ^2 hm-k j-l = °ij ' 

k=0 1=0 
2p Mi i min{AT^j} 

d T y p 5+l3+e + bLPp5 + f3 + e = ~ ^ X P ^ X lS \n ] X ^2 b kl<Pi-k j-l} ■ 

i=p+l j=0 k=0 1=0 

Here [a, bj := [a, 6] n N. This system is easily solved: one begins with i = and 
solves the equations for j running from to Mq. This can then be repeated for 
i = l, etc, until i = p is reached. This provides the functions ipij. Finally, one 
solves the last equation for the remainder term (p p g + p +e , with initial value zero, 
noting that the right hand side of the resulting equation is in 'if So Jr ^ Jr (14 X1 ) , 
and one concludes using Proposition B.l. □. 

PROOF of Proposition B.4: With the notation of the proposition, we set 
fephg = b - bps+e, Cphg = c - Cps+p+e, {p v hg = <P - & P &+f3+t- We use the Lemma 
above to obtain a solution i^phg of the problem 

d T V + bphgP = Cphg , (B.28) 

V Ni 

ip\v = (p = x p ^^2x iS ]n?x(pij{x,v) . (B.29) 
i=0 j=0 

Then we solve 

d T (p' + bip' = c p 5 + /3 +e - bp$ +e (p phg 
with <p'\ T =o = (pps+i3+e- According to Proposition B.l we 

To conclude we set ip = (p p hg + f' which is of the required form, and solves 
(5.24a). ' □ 



B.4 Poly homogeneous solutions of d x (p + bip = c 
Proposition B.6 Let tp be a solution in C^°(U Xl ) of 

d x u> + -<p = c, (B.30) 

x 

and suppose that (B.25) holds with some e e]0,<5[, (3 G R, and with some integer- 
valued sequences (AT/), (iVf). If 

b = o(x) 

(equivalently, boj(0,v,r) = 0), then 

p Mi P Mi 

V = Y.Y. xlS X &3 + x/3+1 E S ^ ^ X + ^+l+/3+6 , (B.31) 
i=0 j=0 i=0 j=0 
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with 

for some integer sequence (Mj). 

Proof: Proposition B.3 shows that for (3 > — 1 the limit 

ifo(-) ■= lim w(x,-) 

exists and is a smooth function on O x [0, T] . When b is a multiple of the identity 
matrix the result is obtained by a straightforward analysis of the formula 

<p(x,-) =<po(-)e-K b (s,-)ds + I* eSW».-)*> c {y,-)dy , (B.32) 

J o 

using the estimates of Lemma B.2. For (3 < —1, and again for b — a multiple 
of the identity matrix — we use instead 

<p(x, •) = <p(x!/2, -)e~ b{s '' )ds + f e% ^'-^cfo, . (B.33) 

Jxi/2 

In the general case, we first note that it follows from Proposition B.3 that there 
exists A6R such that tp G We then write 

^-c = -VC W ; (B.34) 

integrating gives 

Jo 

Inserting this equation in the right-hand-side of (B.34) and integrating again 
one obtains a similar equation with a remainder term falling-off one power of 
8 faster. The result is proved by repeating this procedure a finite number of 
times. □ 
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